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SUMMARY 

We  report  the  results  of  our  theoretical/computational  study  on  the  key  physical  processes 
that  affect  the  propagation  of  high  energy  lasers  in  the  atmosphere.  The  main  objective 
of  this  study  is  to  discuss  the  optimum  laser  wavelength  and  power  for  efficient 
propagation  in  maritime,  desert,  rural  and  urban  atmospheric  environments.  The 
theoretical/numerical  model  used  in  this  study  includes  the  effects  of  aerosol  and 
molecular  scattering,  aerosol  heating  and  vaporization,  thermal  blooming  due  to  aerosol 
and  molecular  absorption,  atmospheric  turbulence,  and  beam  quality.  These  processes  are 
modeled  in  a  fully  three-dimensional  and  time-dependent  manner.  It  is  found  that  aerosol 
particles,  which  consist  of  water,  sea  salt,  organic  matter,  dust,  soot,  biomass  smoke, 
urban  pollutants,  etc.,  are  particularly  important  because  they  result  in  laser  scattering, 
absorption  and  enhanced  thermal  blooming.  In  the  water  vapor  transmission  windows, 
the  total  absorption  coefficient  driving  thermal  blooming  can  be  caused  mainly  by 
aerosols  and  not  water  vapor.  In  certain  maritime  environments  the  deleterious  effects  of 
aerosols  can  be  reduced  by  vaporization.  Aerosol  particles  which  cannot  be  vaporized, 
such  as  those  consisting  of  dust,  soot,  etc.,  can  significantly  increase  thermal  blooming. 
We  show  that  moderate  values  of  the  laser  beam  quality  parameter  have  little  effect  on 
the  propagation  efficiency.  The  laser  power,  averaged  over  dwell  time,  delivered  to  a 
distant  target  as  a  function  of  transmitted  power  is  obtained  for  a  number  of  wavelengths 
and  atmospheric  environments.  The  optimum  wavelength  and  power  are  found  for  each 
atmospheric  environment. 

For  the  particular  propagation  geometry  shown  in  Fig.  t.  our  results  show  that  the 
average  power  on  target  is  strongly  dependent  on  the  atmospheric  environment: 

Maritime:  In  a  maritime  environment,  for  transmitted  power  Pr  <  1 .5  MW.  the 
propagation  efficiency  varies  from  q  *  50%  to  70%.  In  this  transmitted  power  range 
the  1.625  g  m,  and  2. 141  pm  wavelengths  provide  slightly  greater  efficiency  than 
1.045  g  m.  However,  for  PT  >  1 .5  MW.  thermal  blooming  limits  the  power  on  target.  In 
this  high  power  regime,  the  optimum  wavelength  is  1.045  pm  due  to  stronger  absorption 
at  the  other  wavelengths  (see  Fig.  5). 
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Desert:  The  aerosols  in  a  desert  environment  are  composed  mainly  of  dust  particles, 
however,  the  gross  extinction  coefficients  are  similar  to  that  of  a  maritime  environment. 
Hence,  the  power  on  target  for  the  desert  environment  is  very  similar  to  the  maritime 
environment  for  P,  <  1 .5  MW  (see  Fig.  7).  At  higher  transmitted  powers,  however, 
there  is  slightly  less  thermal  blooming  for  1 .045  pm  compared  to  the  maritime 
environment  due  to  the  relatively  lower  aerosol  absorption  in  a  desert  environment. 

Rural:  For  the  rural  environment  2,141  pm  is  the  optimum  wavelength  over  the  entire 
range  of  transmitted  powers.  P,  <  3  MW  .  The  optimum  power  is  found  to  be  Pr—  1.5 

MW  which  results  in  an  average  power  on  target  of  { ~  0.6  MW  and  efficiency  of 

//  *  40% .  For  P,  >1.5  MW.  the  power  on  target  does  not  increase  with  transmitted 
power  because  of  thermal  blooming. 

Urban:  The  optimum  wavelength  for  the  urban  environment  is  found  to  be  2.141  pm 
over  the  entire  range  PT  <  3  MW  (see  Fig.  9).  However,  thermal  blooming  due  to  the 
non-hvgroscopic  aerosols  places  severe  limits  on  the  transmitted  power.  For  example,  the 
maximum  value  of  (fjargel)  ~  0.08  MW  is  obtained  for  Pf  ~  0.2  MW.  giving  an 
efficiency  of  jj  =s  40%  (Fig. 9). 

Air  to  Ground  (Desert,  Urban):  We  have  also  considered  vertical  propagation 
scenarios  from  a  fast  moving  platform  to  a  stationary  target  on  the  ground  in  both  a  desert 
and  urban  environment.  In  these  cases  we  find  that  for  the  desert  environment,  the  large 
beam  slew  negates  thermal  blooming  effects  and  results  in  high  propagation  efficiency. 

//  >  60%  for2.141pm  .  In  the  urban  environment,  however,  the  presence  of  soot 
aerosols  can  still  result  in  significant  thermal  blooming  (see  Fig.  10). 

An  approximate  expression  for  the  Strehl  ratio  of  a  focused  HEL  beam 
propagating  through  a  stagnant  zone  is  derived.  The  propagation  of  a  HEL  beam  in  a 
maritime  atmosphere  with  a  stagnation  zone  is  modeled  using  the  HELCAP  code.  The 
laser  power  delivered  to  the  target  is  calculated  as  a  function  of  slew  rate.  For  the 
parameters  considered,  it  is  found  that  a  stagnation  zone  near  the  laser  source  has  little 
effect  on  the  propagation  efficiency  while  a  stagnation  zone  near  the  target  can 
significantly  reduce  the  power  on  the  target. 
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I.  Introduction 

High  energy  lasers  (HELs)  have  a  number  of  directed  energy  (DE)  applications 
requiring  high-intensity  beams  to  be  propagated  long  distances  under  a  wide  range  of 
atmospheric  conditions.  The  optimum  wavelength  for  efficient  HEL  propagation  depends 
on  the  atmospheric  conditions  and  a  number  of  inter-related  physical  processes  which 
include;  thermal  blooming  due  to  aerosol  and  molecular  absorption  [1].  turbulence  [2], 
aerosol  and  molecular  scattering  [3],  thermal  scattering  due  to  heated  aerosols,  and 
aerosol  heating  and  vaporization  [4-7],  The  relative  importance  of  these  processes 
depends  on  the  parameters  of  the  atmospheric  environment  which  can  vary  significantly 
depending  on  location  and  time. 

Atmospheric  environments  contain  various  types  and  concentrations  of  aerosol 
particles  which  can,  for  HEL  beams,  enhance  thermal  blooming  and  significantly  affect 
the  propagation  efficiency.  In  general,  aerosols  consist  of  hygroscopic  and  non- 
hygroscopic  particles  of  various  sizes  and  chemical  compositions.  Hygroscopic  aerosols 
are  water-soluble  and  vary  in  size  depending  on  the  relative  humidity  [8],  Oceanic 
aerosols  consist  of  sea  salt,  water,  and  organic  material.  Non-hygroscopic  aerosols  are 
composed  of  dust,  soot,  and  other  carbon-based  compounds.  These  aerosols  can  have 
much  larger  absorption  coefficients  than  water- based  aerosols.  While  they  are  normally 
present  in  continental,  rural  and  urban  environments,  dust  aerosol  particles  can  also  be 
present  in  maritime  environments  hundreds  of  miles  from  shore  [9]. 

Aerosols  can  absorb  laser  energy  and.  in  the  case  of  hygroscopic  aerosols,  the 
absorbed  energy  goes  into  both  heating  and  vaporizing  the  aerosol.  Heated  aerosols 
conductively  heat  the  surrounding  air,  resulting  in  an  increase  in  thermal  blooming  of  the 
HEL  beam  [10].  However,  since  aerosol  scattering  and  absorption  coefficients  are 
strongly  dependent  on  the  aerosol  particle  radius,  vaporizing  the  aerosol  can  improve  the 
propagation  efficiency.  Non-hygroscopic  aerosols  (dust,  etc.),  however,  have  large 
scattering  and  absorption  coefficients  and  will  not  vaporize  at  the  intensity  levels 
anticipated  in  DE  applications.  These  aerosols  continually  heat  the  surrounding  air 
leading  to  significant  thermal  blooming. 

Water  vapor  absorption  bands  and  those  of  carbon  dioxide  mainly  determine  the 
atmospheric  transmission  windows  in  the  infrared.  Under  a  range  of  atmospheric 
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conditions  and  laser  wavelengths,  aerosol  absorption  can  exceed  water  vapor  absorption 
and  thus  can  be  the  dominant  process  for  thermal  blooming.  For  example,  in  a  maritime 
environment  at  an  operating  wavelength  of  X  =  1 .045  pm ,  the  water  vapor  absorption 

coefficient  is  -  3x1  O'5  km-1  [11]  while  the  aerosol  absorption  coefficient  is  often  greater 
than  1  O'3  km'1 .  In  other  water  vapor  transmission  windows,  i.e..  1 .625  pm  and  2.141  pm. 
the  water  vapor  and  aerosol  absorption  coefficients  can  be  comparable.  In  addition  to 
enhancing  thermal  blooming,  aerosols  can  also  significantly  contribute  to  the  total  laser 
scattering  coefficient. 

In  this  study,  the  Advanced  Navy  Aerosol  Model  (ANAM)  is  used  to  model  the 
near-surface  maritime  environment  [12].  The  ANAM  aerosol  distribution  is  comprised 
of  various  modes  which  represent  aerosol  particles  of  different  compositions  and  sizes. 
Using  Mie  cross-sections,  we  calculate  the  absorption  and  scattering  coefficients 
associated  with  each  individual  mode.  The  gross  scattering  and  absorption  coefficients 
that  w'e  obtain  are  comparable  with  in  situ  measurements  [13.14], 

In  this  report,  we  analyze  the  relevant  processes  which  limit  HHL  propagation 
efficiency  in  maritime,  desert,  rural,  and  urban  environments.  To  simulate  the  many  inter¬ 
related  processes  affecting  atmospheric  HEL  propagation,  we  use  the  High  Energy  Laser 
Code  for  Atmospheric  Propagation  HELCAP  [15],  developed  at  the  Naval  Research 
Laboratory.  HELCAP  models,  among  others,  the  effects  of  i)  aerosol  and  molecular 
scattering,  ii)  aerosol  heating  and  vaporization,  iii)  thermal  blooming  due  to  both  aerosol 
and  molecular  absorption,  iv)  atmospheric  turbulence,  and  v)  laser  beam  quality.  It  is  the 
first  HEL  propagation  model  w:hich  integrates  all  these  physical  processes  in  a  fully 
three-dimensional,  time-dependent  manner.  In  modeling  the  aerosol  effects,  we  account 
for  the  aerosol  distribution  and  the  various  aerosol  modes  (wrater- based,  dust.  soot.  etc.). 
Furthermore,  since  the  thermal  blooming  process  is  modeled  in  a  fully  time-dependent 
manner,  we  can  simulate  propagation  through  stagnation  zones,  i.e..  locations  at  w'hich 
the  wind/slew  velocity  is  zero  [16]. 

In  Section  II  we  estimate  the  relative  contributions  to  laser  beam  spreading  and 
intensity  loss  in  a  maritime  environment  for  three  wavelengths  lying  within  the  water 
vapor  transmission  window.  The  effects  considered  include:  laser  beam  quality  effects, 
turbulence,  molecular  and  aerosol  thermal  blooming,  and  aerosol  thermal  scattering. 


2 


molecular  and  aerosol  scattering  and  absorption,  In  Section  III.  the  various  aerosol 
models,  e.g.,  the  Navy  Aerosol  Model  (NAM)  [17,18],  the  Advanced  Navy  Aerosol 
Model  (ANAM)  [12],  the  NRL  Aerosol  Analysis  and  Prediction  System  (NAAPS)  [9], 
and  our  method  for  obtaining  aerosol  scattering  and  absorption  coefficients  are  discussed. 
In  Section  IV.  aerosol  heating  and  vaporization  and  their  effect  on  propagation  are 
analyzed.  Thermal  blooming  in  the  presence  of  aerosols  is  analyzed  in  Section  V.  In 
Section  VI.  the  laser  power  delivered  to  a  distant  target  as  a  function  of  transmitted  laser 
power  is  found  for  a  number  of  wavelengths  and  atmospheric  environments,  We  show, 
among  other  things,  that  i)  water  vapor  transmission  windows  are  not  necessarily  the 
determining  factor  for  choosing  the  optimum  HEL  wavelength  ii)  thermal  blooming  due 
to  aerosol  absorption  can  be  the  main  contributor  to  beam  spreading  within  the  water 
vapor  transmission  windows  iii)  non-hygroscopic  aerosols,  because  of  their  large 
absorption  coefficient,  and  the  fact  that  they  can  not  be  vaporized,  are  the  main  sources  of 
aerosol  absorption  and  hence  thermal  blooming  and  iv)  moderate  values  of  beam  quality 
( M 2  <  4 )  have  a  minor  effect  on  the  propagation  efficiency  compared  to  the  effects  of 
turbulence,  thermal  blooming,  and  aerosol  scattering. 

Stagnation  zones  can  be  detrimental  to  HEL  propagation  since,  without  an 
effective  clearing  mechanism  for  the  heated  air,  the  strength  of  the  thermal  lens  grows  in 
time.  In  this  situation,  the  defocusing  process  is  eventually  limited  by  thermal 
conduction  or  buoyancy.  However,  by  the  time  these  processes  become  effective,  the 
beam  may  have  already  been  severely  degraded. 

Thermal  blooming  in  the  presence  of  a  stagnation  zone  was  experimentally 
observed  in  a  laboratory  experiment  using  a  1 OW,  CCh  laser  passed  through  an  absorption 
cell  containing  CO2  gas.  The  cell  was  pivoted  to  simulate  a  wind  profile  containing  a 
stagnation  zone.  This  experiment  was  also  modeled  using  a  code  that  solved  time- 
dependent  thermal  blooming  equations  in  the  isobaric  regime  [16]. 

We  examine  the  effect  of  a  stagnation  zone  on  a  HEL  beam  focused  onto  a  remote 
target  in  a  maritime  atmosphere.  We  first  analyze  the  effect  of  a  stagnation  zone  on  the 
propagat  ion  of  a  focused  beam  and  calculate  an  approximate  expression  for  the  relative 
intensity  (Strehl  ratio)  on  target  as  a  function  of  focusing  geometry  and  stagnation  zone 
position.  We  then  utilize  the  HELCAP  code  [28]  to  examine  the  propagation  of  a 
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megawatt-class  HEL  beam  through  a  maritime  atmosphere  which  contains  a  stationary 
stagnation  zone  along  the  laser  path.  The  laser  target  configuration  used  in  these 
simulations  is  shown  in  Fig.  1.  We  use  the  simulations  to  calculate  the  laser  power 
delivered  to  a  remote  target  for  different  wind  profiles  which  place  the  stagnation  zone  at 
different  locations  along  the  laser  propagation  path. 

II.  Physical  Processes  Affecting  HEL  Propagation 

The  purpose  of  this  section  is  to  obtain  estimates  for  the  relative  importance  of  the 
various  physical  processes  that  lead  to  transverse  spreading  and  loss  of  intensity  of  an 
HEL  beam.  Three  different  laser  wavelengths.^  =  1.045. 1 .625  and  2.141pm .  ail  of 
which  lie  within  water  vapor  transmission  windows,  are  used  for  illustration.  Full  scale 
simulations  of  these  interrelated  processes  will  be  presented  and  discussed  in  Sec.  VI  for 
a  number  of  atmospheric  environments.  In  this  section  we  estimate  these  effects 
individually  in  order  to  better  understand  the  results  of  the  full  scale  simulations. 
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Figure  1:  (a)  Schematic  of  laser  and  target  configuration  used  in  illustration 
and  full  scale  simulations,  (b)  Initial  transverse  intensity  profile  of  the  apertured 
laser  beam  (red  curve)  used  in  the  simulations  with  D  =  80  cm  and  R0  =  50  cm  . 
For  comparison,  the  dashed  curve  denotes  a  Gaussian  beam. 
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Laser  Beam  Quality.  M 2 

4 

Target  Range.  L  [km] 

5 

Wind  Velocity.  VH  [ni/sec] 

5 

Turbulence  Strength.  C2  [m“2  ] 

10_,s 

Water  Vapor  Absorption  Coefficient,  awv  [km-1  ] 

3x1  O'5,  2xl0'\  3xl0'3 

Aerosol  Scattering  Coefficient,  /?A  [km-1] 

1.2xl0-l,7xl0'2.5xl0-2 

Aerosol  Absorption  Coefficient.  ah  [km-1  ] 

2x]  O'1 ,  2x1 0"\  3x1  O'1 

Effective  Aerosol  Absorption  Coefficient,  [km-1  J 

lx  1  O'3 

Table  1:  Laser  and  atmospheric  parameters  used  in  illustration  to 
estimate  and  compare  various  effects 


6 


The  configuration  used  in  the  HEL  propagation  examples  is  shown  in  Fig.  1.  An 
HEL  beam,  with  an  aperture  diameter  of  80cm ,  is  focused  onto  a  target  a  distance  of 
L-  5  km  from  the  source.  For  illustrative  purposes,  we  chose  atmospheric  parameters 
typical  of  a  maritime  environment.  The  formulas  used  in  the  illustration,  however,  can  be 
applied  to  other  atmospheric  conditions.  The  HEL  and  atmospheric  parameters  are  listed 
in  Table  1,  for  three  laser  wavelengths. 

The  average  laser  power  on  the  target  is  determined  by  the  change  in  the  laser 
spot  size  on  target  and  intensity  loss  due  to  the  various  processes.  Some  processes,  e.g. 
turbulence,  are  due  to  small  angle  scattering  events  and  for  our  purposes  are  best 
described  by  a  laser  beam  spreading  angle.  The  HEL  spreading  angle  is  the  ratio  of  the 
change  in  spot  size  to  the  propagation  distance,  i.e..  A©  -  A/? ! L .  Other  processes,  such 
as  molecular  scattering,  result  in  large  angle  scattering  events  and  are  best  described  by 
an  extinction  coefficient. 


i)  Beam  Quality 

It  is  common  practice  to  characterize  the  higher  order  modal  content  of  a  laser 
beam  by  a  beam  quality  parameter  denoted  by  M1 .  The  quantity  A-/2  >  1  is  a  "times 
diffraction-limited"  parameter  which,  for  a  fundamental  Gaussian  beam,  is  unity.  This  is 
one  of  many  measures  of  beam  quality  and  has  a  limited  value  in  determining  the  far  field 
profile.  The  laser  spot  size  on  the  target  due  to  finite  beam  quality,  i.e., 

M 1  =  4.  diffractive  spreading  is 


A  R 


quality 


A 42X 

xRq 


1.1cm,  for  1.045/ym 
<  1.7cm,  for  1.625 /im. 
2.2cm,  for  2.141  jjm 


U) 


ii)  Turbulence 

Temperature  and  density  fluctuations  inherent  in  the  atmosphere  lead  to  random 
fluctuations  in  the  refractive  index.  The  resulting  turbulence  causes  the  laser  beam  to 
transversely  spread  and  wander.  The  size  distribution  of  the  turbulence  is  often  modeled 
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by  a  Kolmogorov  distribution  with  structure  function  parameter  C 2 .  which  characterizes 


the  strength  of  the  turbulence  [2].  The  increase  of  the  laser  beam  spot  size  on  the  target 
due  to  turbulence  is  given  by 


A  R 


turfr 


f  *  \3/5 

clO 


X 


!  1 


2.6  cm,  for  1 .045  pm 
<  2.4cm,  for  1.625pm. 
2.3  cm,  for  2.141pm 


(2) 


Note  that  the  radial  spread.  A Rmrh ,  is  weakly  dependent  on  the  wavelength,  i.e..  is 


proportional  to  A~l's 


iii)  Molecular  scattering 

The  ratio  of  the  laser  intensity  on  a  target  at  range  L  to  that  at  the  source,  due  to 
molecular  scattering,  is 

//souree  =  exp(-/JmI).  (3) 

The  molecular  scattering  coefficient  can  be  written  as 


/?„,  =  n„,  a,„ 

*  m  ffl  iti 


7.5 x  10'4  km'1 ,  for  1. 045 pm 
<  1.3xl0“4  km'1 ,  for  1.625pm. 
4.3x1  O'5  km'1 ,  for  2.141pm 


(4) 


where  nm  is  the  molecular  density  and 

crm  =  ( 8zr / 3 ) (zr («q  -1)/ nm  X  )~  5=  3.3xl0~28  / X(j.un)  is  the  Rayleigh  scattering  cross- 
section.  The  laser  intensity  loss  due  to  molecular  scattering  is  negligible, 

‘•e-’Aarg*  ^source  =  exP {~Pm  L)  *  1  ■  for  all  three  wavelengths. 


iv)  A  erosol  scattering 

The  aerosol  scattering  coefficient  is  fiA  =  J dR  F(R)erscai  (R).  where  F(R)  is  the 

aerosol  particle  radius  distribution  function  and  <Tscat  {R)  is  the  scattering  cross-section 

of  an  aerosol  particle  with  radius  R  .  The  ratio  of  the  laser  intensity  on  a  target  at  range 
L  to  that  at  the  source  is 
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0.5,  for  1 .045  jim 

Aargtf  /  /_  =  exp irfiA  l)  *  <  °-7  >  for  1 . 1 625  pm .  (5) 

0.8,  for  2.141  jim 


Aerosol  scattering  leads  to  significant  loss  of  intensity,  particularly  at  the  shortest 
wavelength. 

v)  Aerosol  thermal  scattering 

Aerosols  absorb  laser  energy  and  heat  the  surrounding  air  through  thermal 
conduction.  The  increase  in  air  temperature  has  a  spatially  fluctuating  component  which 
can  scatter  the  HEL  beam.  The  uniformly  heated  component  of  the  air  temperature 
results  in  thermal  blooming  and  is  discussed  in  subsection  vii).  In  the  geometric  optics 
limit,  multiple,  small-angle  scatterings  result  in  the  spreading  of  the  laser  beam.  The 
increase  in  the  laser  spot  size  on  the  target  due  to  aerosol  thermal  scattering  is  given  by 


^AJ  '  ®AJ  L 


0.16cm,  for  1.045  pm 
n'J2  R*  l)n  *  0.57  cm.  for  1.625  pm  ’ 
1.1cm.  for  2.141pm 

(6) 


*  6.3x1  O'4 


where  ®A  T  is  the  spreading  angle  associated  with  thermal  scattering.  (/)  is  the  average 
laser  intensity  along  the  propagation  path.  Tamh  is  the  ambient  air  temperature.  nA  is  the 
number  density  of  aerosols,  RA  is  the  aerosol  particle  radius,  aD  is  the  bulk  absorption 
coefficient  of  the  aerosols  and  e  is  a  constant  of  order  unity  representing  the  ratio  of 
laser  energy  going  into  vaporization  to  laser  energy  conducted  into  the  air.  In  obtaining 
the  results  in  Eq.  (6)  the  following  values  were  used,  an  =  8.4,  30  and  59  cm-1  for 
1.045  pm,  1,625  pm  and  2.141  pm,  respectively. 

vi)  Thermal  blooming  due  to  water  vapor  absorption 

Molecular  absorption,  particularly  water  vapor  absorption,  heats  the  air  in  the  path 
of  the  HEL  beam  and  results  in  thermal  blooming.  The  molecular  absorption  coefficient 
is  minimized  by  operating  within  the  w'ater  vapor  transmission  window.  The  estimates  in 
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this  and  the  following  subsection  apply  to  whole  beam  thermal  blooming  in  the  steady 
state  isobaric  regime.  A  rough  estimate  for  the  increase  in  spot  size  on  the  target  is 


0.23  cm,  for  1 .045  pm 


«  jlS.Ocm,  for  1.625pm, 
22.5  cm.  for  2.141pm 


(7) 


aUT  is  the  water  vapor  absorption  coefficient, 

Yth  ~  (no  —  \)/p0cpT0  -  7.5x10  '  cm'/J  at  STP,  Vw  is  the  wrind/slewF  velocity,  and 
cp .  p0  and  T0  are  the  specific  heat  at  constant  pressure,  mass  density  and  temperature  of 

air.  respectively.  In  obtaining  Eq.(7)  we  used  the  whole  beam  thermal  blooming 
defocusing  angle  Qm  *  {8nni  R„)L  ®  (nu -\){6p! p„)L! R„  where  Snn  {Sp)  is  the 

variation  across  the  beam  in  the  refractive  index  (air  mass  density)  due  to  thermal 
blooming  as  discussed  in  Sec,  V.  It  should  be  noted  that  in  the  presence  of  wind  or  slew 
the  transverse  intensity  profile  of  the  laser  beam  becomes  highly  asymmetric,  i.e., 
crescent-shaped,  and  the  above  estimate  for  the  spot  size  is  merely  an  indication  of  the 
transverse  scale  associated  with  the  intensity  profile. 

vii)  Aerosol-induced  thermal  blooming 

In  addition  to  the  thermal  scattering  effect  discussed  above,  a  collection  of  heated 
aerosol  particles  can  also  lead  to  enhanced  thermal  blooming  [10].  Aerosol-induced 
thermal  blooming  is  due  to  thermal  conduction  from  the  heated  aerosols  into  the 
surrounding  air.  The  effective  absorption  coefficient  for  aero  sol -induced  thermal 
blooming  is  given  by  a  4  /(]  +  s)  where  aA  is  the  aerosol  absorption  coefficient.  The 
increase  in  the  laser  beam  spot  size  on  the  target  due  to  aerosol-induced  thermal 
blooming  is  approximately 


4.7cm.  fori. 045 um 


(8) 


Aerosol  induced  thermal  blooming  will  be  discussed  in  detail  in  Sec.  V. 

The  contributions  to  the  laser  spot  size  and  loss  in  laser  intensity  on  target  due  to 
the  various  processes  described  above  are  summarized  in  Table  2.  The  laser  and 
atmospheric  parameters  used  in  these  examples  are  listed  in  Table  1  for  three 
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wavelengths  which  lie  within  the  water  vapor  transmission  windows.  The  contribution  to 
the  spot  size  increase  due  to  beam  jitter,  A Rtmr  ~  A0  )Mer  L  ~  1  cm ,  is  the  same  for  the 

three  wavelengths.  Molecular  (Rayleigh)  scattering  is  practically  negligible  in  the  three 
cases. 

Based  on  the  above  illustration  we  find  that:  i)  for  a  laser  wavelength  of  1.045  pm 
the  spread  in  the  beam  spot  size  is  dominated  by  aerosol  induced  thermal  blooming, 
while  the  intensity  on  target  is  reduced  by  almost  50%  as  a  result  of  aerosol  scattering,  ii) 
for  a  laser  wavelength  of  1 .625  pm  thermal  blooming  due  to  water  vapor  absorption  is  the 
dominant  contributor  to  the  spread  in  the  beam  spot  size,  wfhile  the  intensity  on  target  is 
reduced  by  nearly  30%  as  a  result  of  aerosol  scattering,  iii)  for  the  case  of 
2.141  fjm  thermal  blooming  due  to  water  vapor  absorption  is  by  far  the  largest  contributor 
to  the  spread  in  the  beam  spot  size,  while  the  intensity  on  target  is  reduced  by  nearly  20% 
as  a  result  of  aerosol  scattering  and  finally,  iv)  moderate  values  of  the  laser  beam  quality 
factor  M 3 ,  i.e.,  values  less  than  4.  have  little  effect  on  the  propagation  of  HELs  compared 
to  molecular/aerosol  thermal  blooming  effects  or  turbulence. 

Comparing  the  three  wavelengths  considered  in  Tables  I  and  2,  aerosol  scattering 
is  more  important  for  the  shortest  wavelength.  1 .045  pm,  while  water  vapor  induced 

thermal  blooming  is  an  issue  for  the  longest  wavelength.  2.141  pm .  As  far  as  the  loss  in 
intensity  due  to  scattering  is  concerned,  2.141  pm  results  in  the  largest  propagation 

efficiency.  It  should  be  noted,  however,  that  the  results  given  in  Table  2  are  meant  to  be 
illustrative,  and  not  necessarily  typical  of  a  maritime  atmosphere. 

III.  Atmospheric  Aerosols 

As  shown  in  the  previous  Section,  aerosol  scattering  and  absorption  can  play  an 
important  role  in  limiting  the  laser  energy  delivered  to  a  remote  target.  In  typical 
maritime  and  continental  environments,  the  aerosol  scattering  and  absorption  coefficients 
can  be  as  large  as  0.2  km'1  and  0.01km*1 ,  respectively,  even  though  the  average  water 
content  of  aerosols  is  typically  far  less  than  that  of  humid  air.  For  example,  at  a 
temperature  of  30  °C  and  relative  humidity  of  50%,  the  water  vapor  mass  density  is 
pwv  ~  1 .5  x\0~-  g/cm '  while  the  average  mass  density  of  maritime  aerosols  is  typically 


It 


Wavelength,  A  [pm] 

1.045 

1.625 

2.141 

Beam  Quality,  A  R^.  [cm] 

1.1 

1.7 

2.2 

Beam  Jitter,  A Rllller  [cm] 

1 

1 

1 

Turbulence.  A [cm] 

2.6 

2.4 

2.3 

Water  Vapor.  Thermal  Blooming.  A  Rrtt  WV  [cm] 

0.23 

15 

22.5 

Aerosol  Thermal  Scattering.  A  RA  T  [cm] 

0.16 

0.57 

1.1 

Aerosol  Induced  Thermal  Blooming.  A  RTB  A  [cm] 

4.7 

4.7 

7 

Intensity  Ratio  (Molecular  Scattering),  /  lsmrce 

1 

1 

1 

Intensity  Ratio  (Aerosol  Scattering).  !laige,  /  lsnur„ 

0.5 

0.7 

0.8 

Table  2.  Estimates  of  HEL  spreading  and  intensity  loss  due  to  various 
processes  for  three  laser  wavelengths 
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far  less,  <  1 0_)  g/cm? ,  However,  water  molecules  scatter  more  efficiently  in  the  form  of 
aerosols  due  to  the  collective  nature  of  the  scattering. 

Aerosol  particles  occur  over  a  range  of  sizes  and  compositions.  Maritime 
aerosols  consist  of  seawater  droplets  with  radii  in  the  range  0.01  pm  -  10pm  [1 7], 
Continental  aerosols  are  typically  comprised  of  soot  and  non-hygroscopic  dust,  biomass 
smoke,  and  a  variety  of  water-soluble  materials  [9].  There  are  numerous  models  that 
attempt  to  describe  the  size  distribution  and  composition  of  aerosols.  The  Navy  Aerosol 
Analysis  and  Prediction  System  (NAAPS)  is  a  near-operational  predictive  aerosol  model 
that  uses  meteorological  data  from  the  Navy  Operational  Global  Atmospheric  Prediction 
System  (NOGAPS)  to  forecast  aerosol  concentrations  in  real  time.  It  has  extensive 
microphysics  and  chemistry  models  and  includes  dust,  sulfur,  and  smoke  simulations. 

The  Coupled  Ocean-Atmosphere  Mesoscale  Prediction  System  (COAMPS)  is  a  regional 
model  that  works  in  conjunction  with  NAAPS  and  provides  the  vertical  distribution  of 
aerosol  particles  [9],  Presently,  however,  NAAPS  and  COAMPS  have  not  been  applied  to 
specific  near-surface  scenarios  of  interest  for  HEL  applications. 

The  Navy  Aerosol  Model  (NAM)  [17,18]  and  its  successor,  the  Advanced  Navy 
Aerosol  Model  (ANAM)  [12]  are  used  to  model  near-surface  maritime  environments. 
While  ANAM  has  been  benchmarked  in  near-surface  open-ocean  conditions,  it  may  not 
accurately  represent  the  detailed  composition  and  distribution  of  aerosol  particles  in 
regions  where  dust  aerosols  are  expected  to  be  present  [9],  Nevertheless,  ANAM  can 
generate  reasonable  gross  scattering  and  absorption  coefficients  which  are  sufficient  for 
our  purposes  of  simulating  aerosol  induced  thermal  blooming  and  laser  scattering.  In  this 
study,  we  will  use  ANAM  to  generate  the  maritime  aerosol  scattering  and  absorption 
coefficients  used  in  our  simulations. 

The  “Navy  Maritime”  aerosol  model  of  the  MODTRAN  atmospheric 
transmission  code  [19]  uses  NAM,  but  neglects  the  dust  contribution,  i.e.,  mode  0.  Thus, 
it  cannot  be  used  to  accurately  describe  near-shore  maritime  environments.  The  default 
MODTRAN  '‘Maritime”  aerosol  model,  however,  gives  aerosol  absorption  coefficients 
(~KT1  km'1)  that  are  similar  to  ANAM  results  for  polluted  coastal  environments. 
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Mode 

Material 

R.40  M 

Re(rt) 

Im(ff) 

0 

N on-hygroscopic  dust 

0.03 

1.52 

8x1 0'3 

1 

Water  soluable  +  water 

0.03 

1.37 

9.6x1 ,0'5 

2 

Sea  salt  +  water  (“aged"  aerosol) 

0.24 

1.38 

6.9x1 0'5 

3 

Sea  salt  +  water  (new  aerosol) 

2 

1.37 

6.5x1  O'5 

4 

Sea  salt  +  water  (near-surface) 

8 

1.37 

6.5x1 0-5 

Table  3:  Aerosol  material  composition,  mean  radius,  and  refractive  index  of  the  various 
ANAM  aerosol  modes  for  RH  =  80%,  Uto  =  U24  =  5  m/sec,  AMP  =  8.  h  =  5  m. 

X  -  1 .045  pm . 


104 


£  10' 

a 

1 

£ 


u 


10"5 


0-01  0.1  1  10  100 

R  [urn] 

Figure  2:  Aerosol  distribution  function  calculated  according  to  ANAM  3.0.  Dashed 
curves  denote  individual  aerosol  modes  0-4.  Solid  curve  denotes  the  total  aerosol 
distribution  function.  RH  =  80%.  U]0  =  U24  =  5  m/sec.  AMP  =  8,h  =  5m. 

X  =  1 .045  pm . 
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The  ANAM  aerosol  particle  distribution  is  comprised  of  various  modes  which 
represent  aerosol  particles  of  different  compositions  and  sizes.  These  aerosol  modes  will 
absorb  laser  energy  and  vaporize  at  different  rates.  The  aerosol  particle  size  distribution 

4 

function,  F(R)  =  ^  F;  (/J) ,  where  R  is  the  aerosol  particle  radius,  is  represented  as  a 

j  =  o 

superposition  of  five  “modes"  with  each  mode  representing  aerosols  with  a  particular 
physical  composition  and  origin.  The  total  aerosol  density  is  given  by  n A  =  jdR  F(R). 

Mode  0  represents  dust  particles  of  continental  origin,  mode  1  represents  water-soluble 
aerosols,  and  modes  2-4  represent  marine  aerosols  (sea  salt  and  water)  that  result  from 
different  processes.  NAM  contains  only  modes  0-3.  The  physical  properties  of  the 
various  modes  are  summarized  in  Table  3.  Each  mode  is  described  by  a  lognormal 
distribution  over  aerosol  particle  radius  with  a  characteristic  amplitude  and  width.  The 
mean  radius  and  distribution  width  of  the  water-based  modes  (1-4)  are  related  to  the 
ambient  relative  humidity  (RH)  using  the  model  of  Gerber  [20].  The  amplitude  of  modes 
2  and  3  are  related  to  the  24  hour  averaged  wind  speed  (Uj 4)  and  instantaneous  wind 
speed  at  10m  altitude  (C/10),  respectively.  The  amplitude  of  mode  4  is  related  to  the  height 
above  the  sea  surface  ( h )  through  an  empirically  determined  relation.  The  air  mass 
parameter  (AMP)  controls  the  amplitudes  of  modes  0  and  1 .  AMP  is  a  dimensionless 
parameter  varying  between  1  (open  ocean)  and  1 0  (highly  polluted  coastal  area)  that 
qualitatively  characterizes  the  amount  of  dust  or  continental  aerosols  in  the  atmosphere. 
However,  it  is  not  directly  related  to  any  measured  meteorological  parameter  and  can  be 
varied  somewhat  arbitrarily  to  produce  scattering  coefficients  that  agree  with 
measurements.  As  such,  ANAM  has  no  real  predictive  capability  in  regions  w'here  dusty 
aerosols  are  expected  to  play  an  important  role.  Figure  2  plots  the  ANAM  aerosol 
distribution  function  for  the  various  modes  for  RH  =  80%,  £/)0  =  U24  =  5  m/sec,  AMP  = 

8.  and  h  =  5  m. 

Calculation  of  the  aerosol  scattering  and  absorption  coefficients  also  requires  that 
the  complex  refractive  index  of  the  various  aerosol  modes  and  the  complex  cross  section 
of  the  aerosol  droplets,  a  =  crjco,  +  iuahs ,  be  known.  Here,  we  take  the  refractive  index 
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Maritime  Environment  (ANAM) 

Mode 

Density  [cm*3] 

Pa  [km''] 

CL A  [km-'l 

0 

2.6  xlO3 

0.028 

1.6  xlO"3 

I 

6.1  xlO3 

0.038 

4.2  xlO'5 

2 

9.0 

0.032 

3.2  xlO'5 

3 

0.014 

KJ 

X 

O 

1 

u> 

2.4  xlO’5 

4 

0.014 

0.016 

2.2  xlO-4 

Total 

9  xlO3 

0.12 

2  xlO'3 

Table  4:  Aerosol  number  density.  nA ,  scattering  coefficient.  (3 A , 
and  volumetric  absorption  coefficient,  aA  .  associated  with  the 
various  ANAM  aerosol  modes  for  the  same  parameters  as  used  for 
Fig.  2,  i.e„  RH  =  80%.  Uio  =  U24  =  5  m/sec,  AMP  =  8.  h  =  5  m. 

X  =  1.045  pm . 
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of  dust  from  Shettie  and  Fenn  [21],  sea  salt  and  water-soluble  materials  are  taken  from 
from  Volz  [22,  23],  and  pure  water  from  Hale  and  Query  [24].  These  indices  are  also 
tabulated  in  the  technical  documentation  for  NAM  [18].  For  modes  1-4  (hygroscopic 
aerosols)  the  refractive  index  is  also  a  function  of  relative  humidity.  The  values  for  n 
given  in  Table  3  are  calculated  for  RH  =  80  %  [18].  The  complex  cross  sections  are 
calculated  according  to  Mie  theory.  We  use  MODTRAN  to  calculate  the  molecular 
absorption  and  scattering  coefficients  [19]. 

Since  the  aerosol  distribution  can  evolve  with  time,  due  to  vaporization  for 
example,  the  aerosol  absorption  and  scattering  coefficients  are  also  time-dependent  and 
given  by 


f  J* 

<XA(t)= 

/  =  0 

9JiR)Fj(Rtt)dRt 

(9a) 

M)=  {!>,,, 

J=  0 

WJ(tf)F,  (*,/)<«, 

(9b) 

where  the  scattering  and  absorption  cross  sections  are  <rtta(  =  n  Rr  Qwal  ,  and 

,  respectively,  Q  is  the  efficiency  and  j  denotes  mode  number. 

Table  4  lists  the  aerosol  number  density  and  scattering  and  absorption  coefficients 
associated  with  each  aerosol  mode  for  the  same  parameters  used  in  Fig.  2.  For  these 
parameters,  mode  1  (water-soluble  aerosols)  has  the  largest  number  density.  However, 
mode  0  has  the  largest  absorption  coefficient  by  far  due  to  the  large  imaginary  refractive 
index  of  dust-like  aerosols.  Mode  2  has  the  largest  contribution  to  scattering.  The  total 
aerosol  number  density,  scattering  coefficient,  and  absorption  coefficient  associated  with 
the  distribution  of  Fig.  2  are  given  in  Table  2,  nA  =  9x  103  cm'3,  f3A  =  0.12  km"1,  and 
a A  =  t.9x  10"3  km"1 ,  respectively.  The  corresponding  visibility  in  this  example  is 
~3.9  fpA  ~  32km. 


IV.  Aerosol  Heating  and  Vaporization 

Scattering  and  blooming  effects  of  aerosols  can  be  reduced  by  vaporizing  the 
water-based  aerosols.  The  aerosol  scattering  and  absorption  coefficients  are,  in  general. 
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functions  of  the  size  parameter  2  n  RA  /  A .  The  aerosol  absorption  and  scattering 
coefficient  scales  with  aerosol  particle  radius  as, 

aA(0  ~  (10a) 


R6(t),  Rayleigh  limit. 


M)  ~ 


(10b) 


R2{t).  Mielimit, 

where  the  Rayleigh  and  Mie  limits  are  defined  as  2 n  RA  /  A  «  1  and  2 n  RA  /  A  »  1 , 
respectively.  Given  the  strong  dependence  of  orA  and  on  the  aerosol  particle  radius. 


vaporization  can  reduce  both  the  aerosol  absorption  and  scattering  coefficients.  In  the 
following  we  discuss  the  heating  and  vaporization  of  a  single  water- based  aerosol 
droplet,  We  use  these  results  to  model  the  vaporization  of  a  distribution  of  aerosol 
particles  and  the  effect  of  vaporization  on  the  atmospheric  scattering  and  absorption 
coefficients. 


i)  Vaporization  of  an  Aerosol  Droplet 

Heating  and  vaporization  of  a  single  water-based  aerosol  droplet  are  described  by 
the  following  coupled  equations  for  the  aerosol  particle  temperature  and  radius  [6 J. 


dA  TA 
dt 


anI  3  HvopdRr 
+ 


3k 


Pj  ca  ca  «a  Si  p,  ca  R 


t*Ta, 


6  Ra  as  A 


dr 


R 


AT,  , 


(1  la) 


(lib) 


where  A TA  =  TA  -  Tamh  ,  TA  is  the  aerosol  particle  temperature.  Timh  is  the  ambient  air 
temperature.  R is  the  aerosol  particle  radius, 

aD  =  k  R\  Qahs  /(4  x  R]  /  3)  =  3  Qghi  /  4  RA  is  the  bulk  absorption  coefficient  of  the 
aerosol  droplet.  Qiltn(R)  is  the  absorption  efficiency,  as  is  the  evaporation  coefficient 
(sticking  fraction),  pA  is  the  mass  density  of  the  droplet,  k  is  the  thermal  conductivity  of 
air.  cA  is  the  specific  heat  of  the  aerosol  droplet,  H  is  the  enthalpy  of  vaporization. 
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A  =  mv  Da,r  pJexp(~4)/pA  kB  T;mh,  g  =  MvapHvap  /  R  Tamb ,  mt  is  the  weight  of  a  vapor 
molecule.  Doir  =  0.24cnr/sec  is  the  diffusion  coefficient  of  air,  kH  is  the  Boltzman 
constant,  Mvap  is  the  molecular  vapor  mass  (e.g„  Mvap  =18  for  water  vapor). 


R  =  8.3  J  /(K  -  mol)  is  the  universal  gas  constant  and  pa  is  the  constant  of  integration 

(with  units  of  pressure)  in  the  Clausius-Clapeyron  formula,  evaluated  here  for  a  saturated 
(i.e„  100%  RH)  water  vapor  pressure  of  2.34  kPa  at  the  temperature  of  293  K.  The  first 
term  on  the  right  hand  side  of  Eq.  (I  la)  represents  the  absorbed  laser  energy,  the  second 
term  is  due  to  vaporization  and  the  third  is  due  to  thermal  conduction  into  the 
surrounding  air.  The  rate  of  change  of  the  aerosol  particle  radius  is  given  by  Eq.  (1  lb). 
Equations  (1  la)  and  (I  lb)  are  valid  for  £A TA  !TA  «  ]  and 

Tvap  =  !5!n  Ra /af'  »  R]pAcA  Ik.  For  water  at  Tamh  =  293  K ,  cA  =  4.2  J/(g  -  K) , 
k  =  2.5  x  1(T4  W/(cm  -  K)  ,  Hvap  =  2.3  kJ/g  and  we  find  that  £  =  17, 

A  =  2.4 xlO-7  cm'/(K  -  sec) .  For  a  water- based  aerosol  particle  with  R 4  =  1  pm . 

Eqs.(l  la)  and  (lib)  are  valid  for  vaporization  times  rvap  »  0.1  msec.  The  bulk 


absorption  coefficient  for  an  oceanic  aerosol  droplet  is  an  =  8.4,  30  and  59  cm"1  at  the 
wavelengths  A  =  1.045pm,  1.625  pm.  and  2.141pm.  respectively.  Convection  of  the 
aerosols  across  the  laser  beam  due  to  a  wind  or  slew  limits  the  heating  and  vaporization 
time  to  the  local  clearing  time.  This  effect  is  contained  in  the  full-scale  numerical 
simulations  of  Sec.  VI. 

The  aerosol  temperature  increases  due  to  the  absorbed  laser  energy  and  cools  due 
to  vaporization  and  thermal  conduction.  In  the  adiabatic  regime,  where  the  heating  and 
cooling  terms  on  the  right  hand  side  of  Eq.  (11a)  are  balanced,  the  aerosol  temperature  is 
given  by 


AT, 


ah,  i  ft] 

3k-(1  +  e)  ’ 


(12) 


where  e  -  as  Hvap  pA  A  /  k  is  the  ratio  of  the  aerosol  vaporization  energy  to  the  aerosol 
energy  conducted  into  the  air;  i.e„  ratio  of  the  last  terms  in  Eq.(  11a)  [5],  For  water  at  an 
ambient  temperature  of  Tamh  =  293  K  ,  it  is  found  that  e  =  2.2 .  The  adiabatic  regime  is 
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reached  in  a  time  on  the  order  of  the  thermal  conduction  time  given  by 
tm)  =  Paca  Ra  ^(3 /c{1  +  £)).  For  an  aerosol  particle  with RA  =1  pm,  the  thermal 
conduction  time  is  rM  =  20  gsec .  In  the  adiabatic  limit  and  constant  laser  intensity  the 
aerosol  particle  radius  decreases  exponentially  with  time  according  to 

Ra({)  =  RA0  exp(-//rViV),  (13) 

where  ^,t0  is  the  initial  radius  of  the  aerosol  particle  and  t  =  3pAH  (\  +  s~'  )/(aDJ) 
is  the  vaporization  time. 

Figure  3  plots  the  aerosol  particle  temperature  and  radius  as  a  function  of  time  as 

*  •  *  1 

given  by  Eqs.  (11a)  and  (lib).  The  laser  intensity  is  taken  to  be  /  =  2  kW/cm  .  Figure 
(3a)  shows  that  the  temperature  increases  for  a  time  comparable  with  the  thermal 
conduction  time  before  reaching  a  maximum  value  which  is  well-approximated  by  Eq. 
(12).  Figure  (3b)  shows  the  characteristic  exponential  decrease  of  the  aerosol  particle 
radius  in  time.  Note  that  for  these  parameters  the  vaporization  time  is  essentially 
independent  of  initial  radius,  as  predicted  by  Eq.  (13). 


it)  Vaporization  of  a  Distribution  of  Aerosols 

Assuming  that  the  radius  of  each  individual  aerosol  particle  undergoing 
vaporization  evolves  according  to  RA  ( t )  =  RM  h{t ),  as  in  Eq.  (13),  it  can  be  shown  that 


the  aerosol  particle  radius  distribution  function  can  be  written  as 


F(RJ) 


F0(R/h(i)) 

h{t) 


(14) 


where  F0(R)  is  the  initial  distribution  function.  In  the  adiabatic  limit  described  by  Eq. 

(14),  each  mode  is  characterized  by  hft)  -  exp  (-/  /  r  ) ) .  where 

rmp  j  =  3 pA  lHvap  j  (1  +  £-'  )i{aDjl),  for  j  =  1  to  4.  We  assume  that  mode  0  (non- 

hygroscopic  dust)  does  not  vaporize.  In  general.  an  is  weakly  dependent  on  the  particle 
radius.  However,  for  the  purpose  of  obtaining  a  vaporization  time  for  each  mode,  we 
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Figure  3:  Aerosol  temperature  (a)  and  normalized  radius  (b)  versus  time  for  initial 
radii,  RA0  =  0.1, 1,  and  10  pm ,  and  an  I  -  14  kW/cnv  .  Curves  are  almost 
indistinguishable  in  figure  (b). 
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takea,);  -  4 n  Im(n/ )/  X .  which  is  the  absorption  coefficient  in  the  Rayleigh  limit  where 
nt  is  the  refractive  index  for  aerosol  mode  j.  For  a  constant  laser  intensity  of 
1  =  2  kW/cnf  ,  the  vaporization  times  associated  with  each  mode  are  r  ,  =  0.42  sec „ 

and  ^ •■•op. 2  =  ^ iap A  *  0'62  SeC  ■ 

Figure  4  plots  the  total  aerosol  scattering  and  absorption  coefficients  versus  time. 
The  scattering  coefficient  is  seen  to  decrease  by  a  factor  of  5  after  -  1  sec.  The  absorption 
coefficient  is  not  decreased  significantly  due  to  the  presence  of  non-hygroscopic  aerosols 
(mode  0)  which  do  not  vaporize. 


V.  Thermal  Blooming  in  the  Presence  of  Aerosols 


Propagation  of  a  high  energy  laser  beam  in  the  atmosphere  results  in  a  small 
fraction  of  the  laser  energy  being  absorbed  by  both  the  molecular  and  aerosol  constituents 
of  air.  The  absorbed  energy  locally  heats  the  air  and  leads  to  a  decrease  in  the  air  density 
which  modifies  the  refractive  index,  given  by  S nTB  =  (/70  -l)8p/p0,  where  p0  and  dp 

are  the  ambient  and  perturbed  air  mass  densities,  respectively.  The  refractive  index 
variation  leads  to  a  defocusing  or  spreading  of  the  laser  beam  known  as  thermal  blooming 

PI- 

For  an  isobaric  process  the  perturbed  air  temperature  ST  is  related  to  the 
perturbed  density  by  Sp  =  -  (pn  /Tn)ST  and  evolves  in  time  according  to 


c*  Pa 


—  +  yw  v  - 

6t  * 


■\ 

— — —  V2 

cp  J 


&  =  aW!j, 


(15) 


where  k  is  the  thermal  conductivity,  cp  is  the  specific  heat  of  air  at  constant  pressure. 
Vir  is  the  wind  or  slew  velocity  and  I  is  the  time-averaged  laser  intensity.  The  isobaric 
regime  is  valid  for  times  greater  than  the  hydrodynamic  time  R,  / C, ,  where  R,  is  the 
laser  spot  size  and  C\  is  the  acoustic  speed. 


The  rate  of  change  of  laser  energy  density  absorbed  in  air  determines  the  degree 
of  thermal  blooming  and  is  given  by  the  total  absorption  coefficient  aIMal . 
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Figure  4:  Aerosol  absorption  (dashed  curve)  and  scattering  coefficient  (solid 
curve)  for  a  distribution  of  aerosols  versus  time  for  a  constant  laser  intensity, 
1  =  2  kW/cm2 .  The  initial  aerosol  distribution  corresponds  to  that  of  Fig.  2. 
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nnat  '  =  «»r  1  +  4  *  *  »A  RA  A  + 


An 


At 


I *  .amh 


_J  +  ljinA  Pa  *a  A7,  e/?^ 
m..  dt 
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where  the  first  term  on  the  right  hand  side  is  due  to  ambient  water  vapor  absorption,  the 
second  is  due  to  conductive  heating  of  the  air  from  the  heated  aerosols,  the  third  is  due  to 
the  additional  water  vapor  from  the  vaporized  aerosols,  and  the  last  term  is  due  to  the  fact 
that  water  vapor  from  a  vaporized  aerosol  enters  the  air  at  an  elevated  temperature.  The 
third  term  is  small  compared  to  the  first  since  pm.  /  pm  umh  «  1 .  and  the  ratio  of  the 

fourth  to  the  second  term  is 


At  ^  H  RA  d  R  | 

2  mvK  dt 


0.2  at  ATa  !Tamh  «  1. 


Hence,  the  last  two  terms  on  the  right  side  of  Eq.  ( 1 6)  can  be  neglected. 

A  further  simplification  applies  in  the  adiabatic  regime  where  the  temperature 
change  is  proportional  to  the  laser  intensity.  Substituting  Eq.  (12)  for  A T A  into  Eq.  (16) 

results  in  the  following  thermal  blooming  absorption  coefficient  for  a  homogeneous 
aerosol  distribution 


a 


total 


Atr 


"...  ( ' 
i +  £  L  3  j 


(17) 


The  result  in  Eq.(  1 7)  is  important  because  it  shows  that  aerosol  absorption,  modified  by 
vaporization,  contributes  directly  to  the  thermal  blooming  absorption  coefficient. 

In  general,  for  a  non-homogeneous  aerosol  distribution,  the  last  term  in  Eq.(l  7)  must  be 
averaged  over  the  aerosol  distribution  to  give  atolal  =  am  +  a A  /(I  +  £) .  where  we  have 

used  the  definition  of  the  aerosol  absorption  coefficient.  aA  =  jF(/f)£?flAt(/J);r/?  JR 

where  Qghs(R)  =  4  RaJt(R)  /3  is  the  imaginary  part  of  the  scattering  efficiency.  For  the 
multi-mode  aerosol  distribution  of  Fig.  2.  Flq.  (17)  can  be  written  as 
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a 


total 


-  a 


wv 


7  =  0  1  +  Ej 


(18) 
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The  aerosol  contribution  to  the  overall  absorption  coefficient  can  be  much  larger 
than  that  of  molecular  water  vapor.  For  example,  in  the  "water  window"  at  wavelength 
1 ,045pm.  am.  =  3  x  10"'  km"1 ,  while  from  Table  4  the  effective  aerosol  contribution 

can  be  up  to  two  orders  of  magnitude  larger.  When  non-hygroscopic  aerosols  represent  a 
large  fraction  of  the  aerosol  population,  it  is  not  possible  to  significantly  reduce  the 
absorption  coefficient  by  vaporization.  Therefore,  the  optimum  laser  wavelength  for 
reducing  thermal  blooming  should  not  be  primarily  determined  by  the  transmission 
windows  of  molecular  water  vapor,  but  must  also  consider  the  absorption  and  conductive 
air  heating  due  to  aerosols. 

VI.  Simulations  of  HEL  Propagation 

In  this  section  we  present  results  of  full-scale  computer  simulations  of  HEL 
propagation  through  various  atmospheric  environments.  The  propagation  code  used  for 
this  study  is  HELCAP  (High  Energy  Laser  Code  for  Atmospheric  Propagation)  which  is 
a  fully  time-dependent.  3D  code  developed  at  the  Naval  Research  Laboratory'  [15]. 
HELCAP  models  the  propagation  of  continuous  and  pulsed  HELs  through  the 
atmosphere.  Representing  the  laser  electric  field  as 

E  =  A(x,y,z,t) exp[  i(eo0 zlc  -  g)„0]  et/2  +  c.c.,  where  coo  =  2nctX  is  the  laser 
frequency,  eT  is  a  unit  polarization  vector  in  the  x  direction.  A(x,y,z,t)  is  the  complex 

laser  amplitude  and  the  laser  intensity  is  I  =cAA'/%tt.  HELCAP  solves  a  nonlinear 
Schrodinger-like  equation  which  has  the  form 


d_A 

dz 


ic 
2  co.. 


V:A  + 


+  SnTB)  -  ~ {a  +  p) 
c  2 


* 


(19) 


where  a  =  am  +  a A  is  the  total  absorption  coefficient,  /?  =  pm  +  p is  the  total 
scattering  coefficient,  and  dn,  and  SnTB  denote  the  refractive  index  variation  due  to 
atmospheric  turbulence  and  thermal  blooming  respectively.  am (aj  is  the  molecular 
(aerosol)  absorption  coefficient.  ( f$A  )  is  the  molecular  (aerosol)  scattering 
coefficient.  The  quantities^.,  dnTH ,  a  .  and  p  are  space  and  time-dependent  and 
determined  self-consistently  in  the  presence  of  the  effects  discussed  in  the  previous 
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sections,  e.g.,  aerosol  heating  and  vaporization.  The  etYects  of  wind  or  beam  slew,  on  the 
air  and  aerosol  heating  is  contained  in  the  full  scale  simulations  presented  in  this  Section. 

The  terms  denoted  by  represent  other  physical  processes  such  as  group 

J 

velocity  dispersion,  ionization,  relativistic  effects,  nonlinear  Kerr  effects,  and  stimulated 
Raman  scattering.  While  these  processes  do  not  significantly  affect  the  propagation  of  the 
HEL  beams  considered  here,  they  are  important  for  the  propagation  of  ultra-high 
intensity  femtosecond  laser  pulses  [25.  26]. 

In  the  following  examples,  we  consider  the  propagation  of  HEL  beams  in  i) 
maritime,  ii)  desert,  iii)  rural,  and  iv)  urban  atmospheres.  In  these  examples,  the  laser 
wavelengths  are  taken  to  be  1 .045  pm ,  1 .625  pm.  and  2.141  pm.  which  correspond  to 
atmospheric  transmission  windows,  i.e..  minima  in  molecular  (water  vapor)  absorption. 
HELCAP  requires  the  initial  scattering  and  absorption  coefficients  associated  with 
vaporizable  (water-based)  and  non-vaporizable  (e.g..  dust,  soot)  aerosol  constituents  as 
inputs.  For  the  maritime  atmosphere  we  use  AN  AM  to  generate  the  aerosol  scattering  and 
absorption  coefficients.  For  the  urban,  rural,  and  desert  environments  w:e  use 
MODTRAN4  and  the  Air  Force  Geophysics  Laboratory  (AFGL)  model  of  Ref.  [27]  to 
generate  the  aerosol  parameters.  MODTRAN4  yields  the  overall  aerosol  scattering  and 
absorption  coefficients  while  the  AFGL  model  gives  the  physical  compositions  and 
percentages  of  vaporizable  and  non-vaporizable  aerosols. 

The  propagation  configuration  is  shown  in  Fig.l.  The  high  energy  laser  beam  has  an 
initial  field  profile  given  by  A  =  A0  f  (r)g(/)exp(-r2  /  R.I  )  ,  where 

f{r)  =  exp[-(2r /  D )'  ]  ,2  =  20 ,  limits  the  transverse  extent  of  the  beam  to  the  aperture 
diameter.  D.  and  g(t)  is  the  initial  temporal  profile  of  the  beam.  The  transmitted  power 
at  the  source  is  denoted  by  PT  .  The  laser  is  focused  onto  a  remote  target  at  a  range  of  5 
km.  The  target  is  taken  to  be  circular  with  an  area  of  1 00  cm:.  The  propagation  direction 
is  along  the  z-axis  and  a  uniform  transverse  w:ind.  with  velocity  VH  =  5  m/sec  is  directed 
along  the  y-axis.  Atmospheric  turbulence  is  modeled  by  a  Kolmogorov  spectrum  with 
structure  constant  C2  =  10“15  m'23 .  The  pointing  jitter  associated  with  the  laser  beam  is 


26 


taken  to  have  an  angular  spread  of  2  prad  and  a  white  noise  temporal  spectrum.  Since 
thermal  blooming  and  turbulence  can  cause  the  laser  beam  centroid  to  wander,  adaptive 
optics  techniques  are  employed  to  keep  the  laser  beam  centered  on  the  target.  In  the 
simulations,  the  target  is  always  located  such  that  the  peak  laser  fluence  at  5  km  range  is 
at  the  center  of  the  target. 

The  average  power  on  target  is  used  as  a  figure  of  merit  in  the  following 
examples.  It  is  defined  by 


(20) 


where  the  dwell  time  =  1  sec,  and  dxdy  is  the  differential  cross  section  which  is 
integrated  over  the  target  area.  The  total  laser  energy  reaching  the  target  is 


However,  this  laser  energy  is  not  necessarily  absorbed  by  the 


target.  Calculation  of  the  absorbed  laser  energy  requires  additional  information  such  as 
the  target  material  absorption  coefficient,  surface  roughness,  surface  curvature,  etc,  which 
is  not  considered  here. 

i)  Maritime  Environment 

The  maritime  environment  is  characterized  by  a  mixture  of  salt  water  aerosols, 
water  soluble  aerosols  and  dust  aerosols,  as  described  in  Sec.  III.  In  this  example  we  use 
the  aerosol  distribution  shown  in  Fig.  2  to  calculate  the  scanering  and  absorption 
coefficients.  These  coefficients,  as  well  as  the  molecular  absorption  coefficients,  are 
listed  in  Table  1.  For  the  vaporization  calculations,  we  assume  droplet  absorption 
coefficient  of  aI}  =  8  cm"1  ,  30  cm-1 ,  and  59  cm"1  for  the  wavelength  A  =  1.045  pm, 

1 .625  pm  .  and  2.141  pm ,  respectively  .  These  values  are  calculated  assuming  80%  RH 
and  using  the  refractive  index  for  oceanic  aerosols  [18]. 

Figure  5  plots  the  average  power  on  target  versus  the  transmitted  pow'er,  P,  ,  for 
the  three  wavelengths  of  interest.  Our  results  show  that  for  a  maritime  environment,  the 
optimum  wavelength  depends  on  the  transmitted  power.  For  P1  <1.5  MW.  propagation 
is  mostly  affected  by  aerosol  scattering  and  the  average  power  on  target  increases  with 
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Figure  5:  Average  power  on  target,  ( Pf3IiCt } ,  versus  transmitted  power,  Pr .  in  a  maritime 
environment  for  the  wavelengths  A  =  1.045.  1.625.  2.141  pm  .  Initial  beam  profile  has 
R„  =  50  cm .  D  =  80  cm  .  Simulation  geometry  is  shown  in  Fig.  1 .  Aerosol  properties  are 
listed  in  Table  4.  Target  range.  L  =  5  km:  beam  focus  -  5  km:  target  area  =  100  cm"; 
wind  speed.  Vw  =  5  m/sec;  turbulence  strength.  -  10“IS  m'2'3 ;  pointing  jitter  angular 
amplitude  =  2  prad  (white  noise). 
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Pr .  In  this  regime,  the  1 .625  p  m.  and  2.14]  pm  wavelengths  provide  slightly  greater 
power  on  target  than  1 .045  p  m.  This  is  due  to  the  lower  aerosol  scattering  coefficient 
associated  with  the  longer  wavelengths.  For  P,  <1  MW,  the  propagation  efficiency  is 
roughly  50%  for  the  three  wavelengths  considered.  For  example.  Pr  =  1  MW  results  in 
{^target )  ~  0.55  MW  for  X  =  1.625pm  and  X  =  2.141pm.  However,  for  Pt  >  1.5  MW. 

thermal  blooming  becomes  important.  In  this  high  power  regime  the  optimum 
wavelength  is  1 .045  pm  due  to  the  lower  molecular  absorption  coefficient  in  that  water 
vapor  window.  For  PT  =  3  MW,  (P^)  -  I  MW  for  1.045  pm  while  (Pmm)  ~  0.8  MW 

for  1 .625  pm  and  2. 141  pm.  The  power  on  target  decreases  for  larger  values  of 
transmitted  power  not  included  in  the  plot. 

The  effect  of  thermal  blooming  on  the  laser  spot  size  on  target  is  show-n  in  Fig.  6. 
Figure  6  shows  contours  of  the  time  averaged  intensity  in  the  target  plane  for  three  values 
of  transmitted  power  at  X  =  1.625  pm  .  The  time  average  is  performed  over  the  entire 

dwell  time  of  1  sec.  The  laser  spot  size  on  target  is  seen  to  increase  with  increasing  P,  . 
For  relatively  low  power,  PT  =  0.5  MW  ,  the  beam  is  focused  within  the  100  cm2  target 
area  and  the  loss  of  power  on  the  target  is  mainly  due  to  aerosol  scattering.  When 
P,  =  1  MW .  the  laser  beam  extends  slightly  beyond  the  target  area.  For  Pr  =  3  MW ,  the 
beam  cross  section  is  much  larger  than  the  target  area  and  exhibits  a  crescent  shape 
characteristic  of  thermal  blooming  in  the  presence  of  a  wind.  In  these  simulations, 
aerosol  vaporization  effects  increased  the  average  power  on  target  by  ~  1 0% 
for PT  >1.5 MW. 

it)  Desert  Environment 

The  desert  aerosol  environment  is  characterized  by  dry,  dust-like  aerosols  which 
cannot  be  vaporized  at  the  laser  intensities  considered  here.  These  aerosols  absorb  laser 
energy,  heat  the  surrounding  air  and  significantly  contribute  to  thermal  blooming.  The 
aerosol  absorption  and  scattering  coefficients,  generated  using  the  MODTRAN4  desert 
model  with  a  10  km  visibility,  are  shown  in  Table  5  for  the  three  wavelengths  of  interest. 
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Figure  6:  Time-averaged  intensity  in  the  target  plane  for  transmitted  powers,  a)  P 
0.5  MW.  b)  Pj  =  1  MW.  and  c)  Pj  =  3  MW.  and  X  =  1.625  pm  .  Time  average  is 
done  over  1  sec  dwell  time.  Parameters  correspond  to  those  of  Fig,  5. 
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The  molecular  absorption  coefficients  are  taken  to  be  the  same  as  in  the  maritime 
environment. 

Figure  7  plots  the  average  power  on  target  versus  transmitted  power.  PT ,  for  the 
three  wavelengths  of  interest.  The  results  are  qualitatively  similar  to  those  of  the 
maritime  environment,  i.e.,  the  optimum  wavelengths  are  1.625  pm  and  2.141  pm  for 

PT<  2  MW,  and  1.045  pm  for  Pr>  2  MW.  In  a  desert  environment,  PT=  1  MW  results 
in (target)  ~  °-5  MW  for  1-625  pm  and  2.141  pm  wavelengths,  and  (Pwstl)  -  0.35  MW 

for  1.045  pm.  ForPr=3  MW,  (Plslga)  ~  1.3  MW  for  1.045  pm  ,  -  0.8  MW  for  1.625  pm. 
and  -  0.7  MW  for  2.141  pm  wavelengths. 

Hi)  Rural  Environment 

The  rural  aerosol  environment  is  taken  to  be  a  mixture  of  70%  water  soluble 
aerosols  and  30%  dust-like  aerosols  [27].  The  total  aerosol  absorption  and  scattering 
coefficients  for  the  rural  environment  are  shown  in  Table  6.  These  coefficients  are 
generated  using  the  MODTRAN4  rural  model  with  a  10  km  visibility.  At  the  laser 
intensity  levels  considered,  the  dust -like  aerosols  are  not  vaporized  while  the  w:ater 
soluble  aerosols  are  partially  vaporized.  In  the  vaporization  calculations,  the  absorption 
coefficient  of  the  aerosol  droplet  is  taken  to  be  aD  =  1.2  x  10 1  cm"'  ,  1 .0  x  103  cm'1 , 
and  3.5  x  10:  cm'1  for  the  wavelengths  X  =  1.045  pm,  1.625  pm,  and  2.141  pm, 
respectively.  These  values  are  calculated  using  the  refractive  index  for  water  soluble 
aerosols  [27],  assuming  80%  RH. 

Figure  8  plots  the  average  power  on  target  versus  transmitted  power  for  the  rural 
environment.  Because  of  the  large  absorption  coefficient  of  the  water  soluble  and  dust 
aerosols  thermal  blooming  begins  to  be  a  limiting  process  for  P0  >  0.5  MW.  For  the 

rural  environment.  2.141  pm  is  the  optimum  wavelength  over  the  entire  range 
Pr  <  0.3  MW.  The  optimum  power  for  a  wavelength  of  2. 141  pm  is  Pr  ~  1 .5  MW 

which  results  in  ~  0.6  MW.  For  Pr  >  1.5  MW.  the  powder  on  target  is  limited  by 

thermal  blooming. 
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Desert  Environment 

A  [pmj 

[km-1] 

Pa  [km-1] 

1.045 

7xl0"4 

0.17 

1.625 

5x  iO-4 

0.097 

2.141 

6xI0"4 

0.072 

Table  5:  Aerosol  absorption  ( aA )  and  scattering  (ft  A )  coefficients  at 
wavelengths  A  =  1.045, 1.625,  2.141  pm  for  a  model  desert  environment 
generated  using  MODTRAN4.  10  km  visibility. 


Figure  7:  Average  power  on  target.  ) .  versus  transmitted  power,  Pr ,  in  a 
desert  environment  for  the  wavelengths  A  =  1.045, 1.625.  2.141  pm  .  Initial  beam 
profile  has  =  50  cm .  D  -  80  cm  .  Simulation  geometry  is  shown  in  Fig.  1 . 
Aerosol  properties  are  listed  in  Table  5.  Target  range.  L  =  5  km;  beam  focus  =  5  km; 
target  area  =  1 00  cm2;  wind  speed.  Vw  -  5  m/sec:  turbulence  strength. 

C';t  =  10“h  nT" 3 ;  pointing  jitter  angular  amplitude  =  2  prad  (white  noise). 
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iv)  Urban  Environment 

The  aerosol  absorption  and  scattering  coefficients  are  generated  using  the 
MODTRAN4  urban  aerosol  model  with  a  10  km  visibility.  Table  7  lists  the  calculated 
scattering  and  absorption  coefficients.  The  aerosol  absorption  in  an  urban  environment  is 

the  largest  of  the  four  environments  considered.  Urban  aerosols  are  modeled  as  a  mixture 
of  80%  rural  aerosols  and  20%  soot  aerosols  [27].  Soot  aerosols,  which  cannot  be 
vaporized,  represent  the  dominant  contribution  to  aerosol  absorption.  Hence  they  heat  the 
air  and  cause  significant  thermal  blooming  of  the  laser. 

Figure  9  plots  the  average  power  on  target  versus  the  transmitted  power.  For  the 
urban  environment.  2.141  pm  is  the  optimum  wavelength  over  the  entire  range 

Pr  <  3  MW.  The  optimum  transmitted  power  for  the  2.141  pm  wavelength  is  Pr  ~  0.3 

MW.  which  results  in  a  power  on  target  of  only  (F^,)  -  0.09  MW. 

v)  Air-to-Ground  Propagation  in  Desert  and  Urban  Environments 

Finally,  we  consider  the  propagation  of  a  H EL  beam  from  a  high  altitude  fast- 
moving  airborne  platform  or  plane,  to  a  stationary  target  on  the  ground.  The  plane  is 
located  at  z  =  0  and  the  target  is  at  z  =  L  =  5  km.  The  laser  beam  has  an  effective  slew 
velocity  of  Vw  =  V0[l  -  (z/Z,)],  where  V0  =  100  m/sec  is  the  plane  velocity.  The 

background  wind  velocity  is  zero.  This  configuration  produces  a  stagnation  zone  near  the 
ground,  i.e.,  target.  The  turbulence  strength,  scattering  coefficients,  and  absorption 
coefficients  are  assumed  to  vary  with  atmospheric  density  according  to,  for  example. 

Cl  -  C„2  g  exp[(z  -  L)  /  Lm  ]  where  C2n  g  =  10'1S  m'33  is  the  turbulence  strength  at 
ground  level  and  Lam  =  8  km  is  the  characteristic  height  scale  for  the  atmospheric 

density.  All  other  parameters  are  the  same  as  in  the  previous  examples.  The  range  to  the 
target  is  assumed  to  be  constant. 

Figure  10  plots  the  average  power  on  target  versus  the  transmitted  power  for  a 
desert  and  urban  environment.  It  is  seen  that  in  the  desert  environment  the  large  slew' 
effectively  reduces  thermal  blooming  and  results  in  propagation  efficiencies  of  greater 
than  rj  %  60%  for  the  case  of  the  2.141  pm  wavelength  laser.  In  the  urban  environment. 
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Rural  Environment 

X  [pmj 

[km'1] 

A  I*™'1! 

1.045 

0.016 

0.15 

1.625 

0.012 

0.076 

2.14! 

0.006 

0.053 

Table  6:  Aerosol  absorption  ( aA )  and  scattering  ( (3A )  coefficients  at 
wavelengths  X  -  1.045, 1.625,  2.141  pm  for  a  model  rural  environment 
generated  using  MODTRAN4.  10  km  visibility. 


Figure  8:  Average  power  on  target.  ) ,  versus  transmitted  power.  P,  .  in  a  rural 
environment  for  the  wavelengths  X  =  1 .045.  1 .625.  2.141  pm  .  Initial  beam  profile 
has  Rn  =  50  cm .  D  ~  80  cm  .  Simulation  geometry  is  shown  in  Fig.  1.  Aerosol 
properties  are  listed  in  Table  6.  Target  range.  L  =  5  km:  beam  focus  =  5  km;  target 
area  =  100  cm2;  wind  speed.  Vw  =  5  m/sec;  turbulence  strength.  C2  =  1 0” 15  m'23 ; 
pointing  jitter  angular  amplitude  =  2  prad  (white  noise). 
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Urban  Environment 

k  [pm] 

aA  [km'1] 

Pa  [km'1] 

1.045 

0.05 

0.13 

1.625 

0.036 

0.065 

2.141 

0.028 

0.044 

Table  7:  Aerosol  absorption  (aA )  and  scattering  ( )  coefficients  at 
wavelengths  k  -  1.045, 1.625,  2.141  pm  for  a  mode!  urban  environment 
generated  using  MODTRAN4.  10  km  visibility. 


Figure  9:  Average  power  on  target.  {PW9n ) ,  versus  transmitted  power.  P,  ,  in  an 
urban  environment  for  the  wavelengths  k  =  1.045. 1.625,  2.141  pm  .  Initial  beam 
profile  has  Rg  =  50  cm ,  D  =  80  cm  .  Simulation  geometry  is  shown  in  Fig.  3 . 

Aerosol  properties  are  listed  in  Table  7.  Target  range,  L  =  5  km;  beam  focus  =  5  km; 
target  area  =  100  cm2;  wind  speed,  Vw  =  5  m/sec;  turbulence  strength. 

Cl  =10'"  in'" 3 ;  pointing  jitter  angular  amplitude  =  2  prad  (white  noise). 
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thermal  blooming  is  still  a  limiting  factor  due  to  the  large  absorption  of  soot  aerosols. 
However,  the  propagation  efficiency  in  this  vertical  propagation  example  is  much  greater 
than  for  the  urban  horizontal  propagation  example  of  Fig.  9.  For  the  optimum  wavelength 
of  2.141pm  ,  the  propagation  efficiency  at  PT  =  1  MW  is  //  =  45%. 
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Transmitted  Power,  PT  [MW] 


Figure  1 0:  Average  power  on  target.  ) ,  versus  transmitted  power,  Pr ,  in  a  a) 

desert  environment  and  b)  urban  environment  for  the  wavelengths  A  =  1 .045. 1 .625, 
and  2.141  pm  for  vertical  air-to-ground  propagation,  Initial  beam  profile  has 

=  50  cm ,  D  =  80  cm  .  Simulation  geometry  is  shown  in  Fig.  1 .  Aerosol  properties 
are  listed  in  Tables  5  and  7.  Target  range,  L  =  5  km;  beam  focus  =  5  km;  target  area  = 
100  cm2;  wind  speed.  Vw  -  V0[l  -  (z/ £)].  Vw  =  100  m/sec;  turbulence  strength, 

=Cl.g  exp[(z-I)/I0,m],  C„2g  =10'"  m'2'\  Lalm  =  8  km ;  pointing  jitter  angular 
amplitude  =  2  prad  (white  noise). 
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VII.  HEL  Propagation  Through  a  Stagnation  Zone 

High  energy  laser  beams  propagating  through  the  atmosphere  can  be  severely 
defocused  or  deflected  by  thermal  blooming  [1  ].The  thermal  blooming  process  is  driven 
by  a  small  fraction  of  the  laser  energy  that  is  absorbed  by  the  molecular  and  aerosol 
constituents  of  air  [10,28],  The  absorbed  energy  locally  heats  the  air  and  leads  to  a 
decrease  in  the  air  density  which  modifies  the  refractive  index.  The  refractive  index 
variation  leads  to  a  defocusing  or  deflection  of  the  laser  beam.  In  the  presence  of  a 
transverse  wind,  the  region  of  heated  air  is  convected  out  of  the  beam  path  and  a  steady- 
state  situation  is  realized  [  1  ].  In  general,  however,  the  intensity  of  a  beam  undergoing 
thermal  blooming  is  a  function  of  both  time  and  spatial  position,  particularly  in  a 
stagnation  zone,  where  the  effective  wind  velocity  is  zero. 

Stagnation  zones  are  particularly  detrimental  to  HEL  propagation  since,  without 
an  effective  clearing  mechanism  for  the  heated  air.  the  strength  of  the  thermal  lens  grows 
in  time.  In  this  situation,  the  defocusing  process  is  eventually  limited  by  thermal 
conduction  or  buoyancy.  However,  by  the  time  these  processes  become  effective,  the 
beam  may  have  already  been  severely  degraded. 

Thermal  blooming  in  the  presence  of  a  stagnation  zone  was  experimentally 
observed  in  a  laboratory  experiment  using  a  1 OW.  CO;  laser  passed  through  an  absorption 
cell  containing  CCL  gas.  The  cell  was  pivoted  to  simulate  a  wind  profile  containing  a 
stagnation  zone.  This  experiment  was  also  modeled  using  a  code  that  solved  time- 
dependent  thermal  blooming  equations  in  the  isobaric  regime  [16]. 

In  this  Section  we  examine  the  effect  of  a  stagnation  zone  on  a  HEL  beam 
focused  onto  a  remote  target  in  a  maritime  atmosphere.  We  first  analyze  the  effect  of  a 
stagnation  zone  on  the  propagation  of  a  focused  beam  and  calculate  an  approximate 
expression  for  the  relative  intensity  (Strehl  ratio)  on  target  as  a  function  of  focusing 
geometry  and  stagnation  zone  position.  We  then  utilize  the  HELCAP  code  [28]  to 
examine  the  propagation  of  a  megawatt-class  HEL  beam  through  a  maritime  atmosphere 
which  contains  a  stationary'  stagnation  zone  along  the  laser  path.  The  laser  target 
configuration  used  in  these  simulations  is  show:n  in  Fig.  1.  Our  simulation  model  includes 
time-dependent  thermal  blooming  as  w:ell  as  aerosol  and  molecular  absorption  [  1 0,28], 
turbulence  [2],  aerosol  and  molecular  scattering  [3],  and  aerosol  heating  and  vaporization 
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[4-6],  We  use  the  simulations  to  calculate  the  laser  power  delivered  to  a  remote  target  for 
different  wind  profiles  which  place  the  stagnation  zone  at  different  locations  along  the 
laser  propagation  path. 


i)  Analysis 

The  effect  of  a  stagnation  zone  on  HEL.  propagation  in  the  atmosphere  is  highly 
dependent  on  laser  and  atmospheric  parameters  as  well  as  on  the  propagation  geometry. 
In  general,  the  problem  requires  a  numerical  simulation  to  determine  the  amount  of 
power  reaching  a  remote  target.  Before  presenting  results  of  full-scale  propagation 
simulations,  however,  it  is  useful  to  analyze  the  following  simpler  problem. 

Consider  a  focused  laser  beam  of  wavelength,  A ,  propagating  from  z  =  0  to  a 
remote  target  at  a  range  z  =  L,  where  z  denotes  the  propagation  coordinate.  The  focal 


length  of  the  laser  beam.  Lf ,  is  such  that  Lf  >  L.  The  wind  speed,  V(z)  is  only  a 
function  of  z  and  has  a  stagnation  point  defined  by  V(zQ)  =  0,  where  0  <  z0  <  L  .  The 

atmosphere  is  characterized  by  an  absorption  coefficient,  or ,  and  a  scattering  coefficient. 
P .  Turbulence  is  neglected  in  this  analysis.  This  configuration  allows  an  approximate 
expression  for  the  Strehl  ratio  (relative  intensity)  at  the  target  to  be  calculated  analytically 
as  a  function  of  both  the  location  of  the  absorbing  region  and  the  laser  focal  length. 

The  propagation  of  the  laser  beam  is  described  by  the  paraxial  wave  equation 


dA 

dz 


—  V*  +  /  —SnTB  -  ^ 
2k„  1  c  ™  2 


A , 


(21) 


where  the  laser  electric  field  is  E  =  A  exp[  i{coa  zl  c  —  co0 1)]  e(  /  2  +  c.c.,  (Ot,  -  2nd  A . 
et  is  a  unit  polarization  vector  in  the  x  direction,  k0  =  n0  (o0 ! c ,  n0  is  the  ambient 


refractive  index,  V1  is  the  gradient  operator  in  the  transverse  direction,  c  is  the  speed  of 
light  in  vacuum,  Snm  is  the  change  in  refractive  index  induced  by  the  laser,  and  c.c. 
denotes  the  complex  conjugate.  In  the  absence  of  a  transverse  wind  or  slew,  there  is  axial 
symmetry'  about  the  z-axis  and  all  spatially  varying  quantities  are  functions  of  (r,z,/) 
only,  where  r  is  the  radial  coordinate  and  t  is  time.  Writing 

A(r,z,t)  =  fi(r,z./)exp[/0(r,z,/)],  where  B ,  6.  and  8  nTB  are  assumed  to  be  real,  leads 
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to  the  following  equations  for  the  laser  intensity  /  (~  B~ )  and  the  transverse 
wavenumber  k  =  V±  6  : 


8J_ 

8z 


-fv±.(k  J)-pi. 

k„ 


(22) 


gk, 

8z 


=  |i-Vx|-l(2/Vl/  -  (V,/)2)!  -  -pk1.V1k1  +  (23) 


On  the  right  side  of  Eq.  (23).  the  first  term  is  of  order  ~  A  i  R  '  where  R  is  the  laser  spot 
size.  The  second  term  is  of  order  -  R/(L)  A)  for  a  focused  beam  and  the  third  term  is  of 
order  ~  5  nm  /(A  R) .  We  consider  parameters  such  that  R  ~50  cm ,  A  -  1  pm ,  and 
L{  —  5  km  .  In  this  parameter  regime,  the  second  term  is  much  larger  than  the  first  term. 


The  third  term  is  much  larger  than  the  second  term  for  dnm  >  10  s .  w'hich  is  typically 
encountered  within  a  stagnation  zone.  Hence,  for  the  parameter  regime  of  interest,  we 
neglect  the  first  and  second  terms  on  the  right  side  of  Eq.  (23)  and  integrate  to  solve  for 
kx,  i.e.f 


co  "V 

kL(/\z.r)*  ki0(r)  +  —  \V±&irg(r,z',t)dz', 

c  J 


(24) 


where  kx0(r)  is  the  initial  transverse  wavenumber  at  z  =  0.  Equation  (24)  represents 
k  in  the  geometric  optics  limit. 

Equation  (22)  can  be  integrated  to  yield 

I(r,z,t)=  exp]  — —  f- — - — Vx  •  [kx(r,z’,/)/(/\ z1. /)]<£'  ,  (25) 

[  w0oJ/(r,z,/)  j 

w'hich  is  an  integral  equation  for  I(r,z,l)  that  can  be  solved  approximately  by  replacing 

I(r.zj)  in  the  right  side  by  /0(r,z) .  i.e..  the  unperturbed  intensity  for  a  homogeneous 

medium  ( S  nIH  =  0 ).  The  unperturbed  intensity  satisfies 


/0(r.z)  =  I(t\Q)e~P:  expj — —  J  1  Vx  -  [k10(r)/0(f,z')]<fe' 

1  o>0  ih(r,z)  J 

Substituting  Eq.  (24)  into  Eq.  (25)  yields  the  relative  intensity  [29] 


(26) 
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The  refractive  index  change  due  to  thermal  blooming  is  given  by 
Sn Ta  =  («„  -\)8p! Pq  .  where  p0  is  the  mass  density,  and  5p  is  the  perturbed  mass 
density.  In  the  isobaric  regime,  the  mass  density  change  is  given  by  [1,28] 


—  +  V-V  - 


— V2  8p  -  -aTH(r,t)I , 


C  T 

'-'P  l0 


(28) 


dr 


^  P  Po  ) 


where  at  is  the  thermal  conductivity,  V  is  the  effective  wind  velocity,  C  is  the  specific 
heat  at  constant  pressure,  T0  is  the  ambient  temperature,  and  I  is  the  time-averaged  laser 
intensity.  In  general,  the  absorption  coefficient.  aw ,  contains  both  aerosol  and  molecular 
contributions  and  can  be  time  dependent.  For  the  simplified  analysis  of  this  Section, 
however,  we  take  it  to  be  constant  in  time.  The  isobaric  approximation  is  valid  for  times 
greater  than  the  hydrodynamic  time  R/C%.  where  C(  is  the  acoustic  speed.  We  consider 

time  scales  less  than  the  thermal  conduction  time,  C  p0R2  Ik ,  and  neglect  the  thermal 
conduction  term  in  Eq.  (28). 

For  the  purpose  of  obtaining  an  approximate  closed  form  expression  for  Eq.  (27), 
we  solve  Eq.  (28)  approximately  by  writing  the  operator  V-V  a  1  /  rt.  (z)  where 

rt.(z)  =  /?(z)/|  V(z)\  is  the  clearing  time.  Using  this  approximation,  we  can  integrate  Eq. 
(28)  to  obtain 


8nm(r,z,t)  »  -ynaI0(r,z)Te{z){ 1  -  exp[- / / rf (z)]} , 


(29) 


where  ym  =  (n0  - \)l(CpT0p0) . 


The  unperturbed  intensity  is  taken  to  have  a  Gaussian  transverse  profile 


(30) 
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where  denotes  the  intensity  on  axis  at  z  =  0 .  The  spot  size  R  is  a  function  of  z  and 
the  initial  spot  size  is  denoted  by  R0  =  R(z  =  0) .  Substituting  Eqs.  (29)  and  (30)  into  Eq. 
(27),  resuits  in  an  expression  for  the  on-axis  Strehl  ratio  at  z  =  L ,  i.e.. 


/(r  =  0.  L,t) 
I0(r  =  Q.L) 


exp 


i.  * 

-8/m  RlYwa\dz'\dz"G{z",t) 


(31) 


where  G(z,t)  =  rr(z)[l  -  exp(-r/rt  (z))]exp(-/?z)/  R*(z).  To  obtain  an  approximate 

analytic  expression  for  the  right  side  of  Eq.  (3 1 ),  we  note  that  for  a  weakly  focused  beam, 
the  integrand  G  is  sharply  peaked  and  grows  linearly  with  time  in  the  vicinity  of  the 
stagnation  zone.  i.e..  G(z0.t)  =  /exp(-/?z0)/  /?4(z0).  Hence,  we  can  approximate  the 
integration  over  z"  in  Eq.  (31 )  by 


\dz"  G(z”,t)  *  G(z0,OAz©(z’-z0).  (32) 

o 

where  A z  is  the  characteristic  width  of  the  stagnation  region  and  0  is  the  Heaviside  step 
function. 

Using  Eqs.  (32)  in  Eq.  (31)  results  in 


/  (r  =  0  ■/../) 

A>0*  =  <U) 


*  exp 


8y/wa/ oaL2t  f  AzY 


R<\ 


1  _  £oA*o4  exp(-/?z0) 

J 


(33) 


The  width  of  the  stagnation  region,  Az,  is  in  general  a  function  of  time  and  also 
the  spatial  variation  of  the  wind  profile.  We  can  define  the  stagnation  zone  as  the  region 
where  the  condition  1 1  rt.(z)  «  1  is  satisfied.  In  this  case,  the  width  of  the  stagnation 

zone  can  an  be  defined  according  to  t  /  r(,(z0  +  Az)  =  1 . 

We  now  consider  specific  analytic  expressions  for  the  wind  profile  V (z)  and  laser 
spot  size  R(z) .  For  a  slewed  laser  beam,  the  wind  profile  can  be  written  as 


nz)  =  n-d-z/z.).  (34) 

where  Vwmd  is  the  ambient  wind  velocity,  and  the  laser  slew  rate  .  0  is  related  to  z0 
according  to  0  =  Vwmd  /  z0 .  Note  that  for  this  model,  a  stagnation  zone  at  the  transmitter 
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( z  =  0 )  implies  an  infinite  slew  rate.  Assuming  that  the  laser  spot  size  does  not  vary 
appreciably  within  the  stagnation  zone. 


At  *  z0R(z0)/(Vwmdt) 


(35) 


for  the  wind  profile  of  Eq.  (34). 

For  the  linearly  varying  wind  profile  described  above.  Eq.  (33)  can  be  written  as 


(36) 


where 


(37) 


describes  how  the  location  of  the  stagnation  zone  affects  the  relative  intensity  on  target. 
The  larger  the  value  of  f{zQ),  the  lower  the  intensity  on  the  target. 

We  can  derive  an  approximate  functional  form  for  R(z)  to  be  used  in  Eq.  (37). 
Substituting  Eq.  (30)  into  Eq.  (26),  the  laser  spot  size  can  be  shown  to  satisfy 
R(z)  =  R0[]  -  z(\  -  rj)/ L]  for  z  «  Lf  ,  where  the  parameter  rj  s  R(L)/ R0,  i.e.,  the 

ratio  of  the  spot  size  on  target  to  the  initial  spot  size.  In  terms  of  the  focal  length. 

t}  =  1  -  (£/ Lf).  Phenomenologically,  rj  can  also  parameterize  focusing  limitations  due 

to  beam  quality  or  atmospheric  turbulence. 

Figure  1 1  plots  the  distortion  function  /  versus  z0  for  various  values  oft;.  We 
assume  a  normalized  scattering  coefficient  of  (3  L  =  0.25 .  For  a  beam  that  is  focused 
near  the  target  plane,  the  distortion  at  the  target  plane  is  largest  when  the  stagnation  zone 
is  close  to  the  target.  In  contrast,  when  the  beam  is  collimated  ( rj  ~  1 ),  the  distortion  at 
the  target  plane  is  greater  when  the  stagnation  zone  is  closer  to  the  laser  source.  Note  that 
the  case  z0  =  0,  which  corresponds  to  an  infinite  slew  rate  given  the  assumed  wind 
profile  of  our  model,  is  not  physically  realizable. 
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Figure  1 1 :  Distortion  function  /  as  defined  by  Eq.  ( 1 7)  plotted  versus  distance  to 
stagnation  zone  zn  for  various  values  of  the  focusing  parameter  q  =  R(L)/  R0  where 
R(L)  is  the  unperturbed  laser  spot  size  at  the  target  plane  ( z  =  L ),  Ra  is  the  initial  spot 
size  at  z  =  0 ,  and  p  L  =  0.25 . 

it)  Numerical  Simulation 

For  the  numerical  simulations  presented  here,  we  use  the  HELCAP  code  and 
consider  a  slewed  HEL  beam  focused  at  a  range  of  L  =  5  km  in  the  presence  of  a  uniform 
transverse  wind.  The  z-axis  defines  the  direction  of  propagation  and  the  wind  and  slew- 
velocities  are  taken  to  be  along  the  y-axis.  The  simulation  geometry  is  shown  in  Fig.  1 . 
The  high  energy  laser  beam  has  an  initial  field  profile  given  by 
A  =  T0/(r)g(/)exp(-r2  /  ) ,  where  f(r)  =  exp[-(2r/  D)'  ]  ,£  =  20,  limits  the 

transverse  extent  of  the  beam  to  the  aperture  diameter,  D.  The  function  g(/)  is  the  initial 
temporal  profile  of  the  beam  which  contains  a  pointing  jitter  characterized  by  an  angular 
spread  of  2  prad  and  a  white  noise  temporal  spectrum.  The  transmitted  power  at  the 
source  is  denoted  by  P,  . 
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The  laser  is  propagated  through  a  maritime  environment  in  which  the  aerosol 
distribution  is  modeled  using  the  Advanced  Navy  Aerosol  Model  (ANAM)  [12]  and  the 
molecular  absorption  coefficient  is  calculated  using  MODTRAN4  [19].  The  details  of  the 
atmospheric  model  and  the  parameters  used  to  characterize  a  maritime  environment  are 
discussed  in  previous  chapters  of  this  report.  The  ambient  absorption  and  scattering 
coefficients  for  propagation  at  A  =  2.141  /aware  am  =  6x10'*  km'1  and 

/?  =  0.05  km"1 ,  respectively,  In  the  simulations  these  quantities  vary  in  space  and  time 
due  to  aerosol  vaporization. 

The  effective  wind  speed  along  the  propagation  path  is  V(z)  =  Vwmd  -  9z.  where 
z  is  the  distance  from  the  laser  source.  Vwnd  is  the  ambient  wind  speed,  and  9  is  the  laser 

beam  slew  rate.  For  positive  9 ,  i.e.,  slewing  in  the  direction  of  the  wind,  a  stagnation 
zone  is  created  at  location  z0  =  Vwnd  19 .  In  what  follows,  we  investigate  the  effect  of 

this  stagnation  zone  on  the  propagation  of  a  HEL  beam  for  various  slew  rates  which 
change  the  distance  from  the  transmitter  to  the  stagnation  zone. 


Slew  [mrad/sec] 


Figure  12:  Average  power  on  target,  {  Flargel) .  versus  slew  rate  for  X  =  2.141pm, 

PT  =  1 .5  MW  ,  VwM  =  5  m/sec,  tdweU  =  1  sec ,  and  L  =  5  km .  Points  (a),  (b).  (c),  and  (d) 
are  labeled  to  correspond  with  fluence  distributions  shown  in  Fig.  4. 
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iii)  Results 

We  consider  a  laser  of  wavelength  X  =  2.141pm  and  power  PT  =1.5  MW 
propagating  through  a  maritime  atmosphere  with  Vwmd  =  5  m/sec.  The  power  delivered 

to  a  target  of  circular  area  1 00  cirr  at  a  range  of  5  km  is  used  as  a  figure  of  merit.  In  the 
simulations  that  follow,  the  target  is  always  located  such  that  the  peak  laser  fluence  at  5 
km  range  is  at  the  center  of  the  target.  The  average  powrer  reaching  the  target  is  defined 
by 

1  Tj"r 

= - \dr  jdxctyI(x,y,z  =  L,T),  (38) 

Zjwell  o 

where  the  dwell  time  =  1  sec ,  and  dxdy  is  the  differential  cross  section  which  is 
integrated  over  the  target  area.  The  total  laser  energy  reaching  the  target  is 
E  =(P  \t 

-  target  \ 4  target  /  4  ttwett  * 

Figure  12  plots  (/>lar^el )  versus  slew  rate,  9 .  for  the  propagation  configuration  of 

Fig.  1 .  A  negative  slew  denotes  that  the  slew  direction  is  opposite  to  the  wind  direction, 
hi  this  case,  there  is  no  stagnation  zone  between  the  laser  and  the  target.  A  positive  slew' 
denotes  slew  in  the  direction  of  the  wind.  In  this  case,  it  is  possible  to  have  a  stagnation 
zone  somewhere  along  the  propagation  path.  For  example,  a  slew  rate  of  1  mrad/sec 
places  the  stagnation  zone  at  the  target,  while  a  slew1  rate  of  5  mrad/sec  places  the 
stagnation  zone  1  km  from  the  laser  source.  Figure  12  shows  that  in  the  absence  of  slew 

(i.e..  uniform  wind),  the  power  on  target  is  (Px^a)  ~  0.7  MW  and  the  propagation 

efficiency  is  -  50%.  ( Paatet )  increases  relative  to  this  value  when  the  slew'  is  opposite  to 

the  wind  direction  because  the  effectively  larger  wind  speed  along  the  entire  propagation 
path  mitigates  thermal  blooming.  For  a  slew^  rate  9  = -5  mrad/sec .  for  example. 

(^targei )  *  1  ■  1  MW .  which  corresponds  to  a  propagation  efficiency  of  >  70%. 

When  the  slew  is  increased  in  the  direction  of  the  wind,  the  power  on  target 
decreases  sharply  for  6  <  l  mrad  / sec  and  reaches  a  minimum  when  9  -  1  mrad  /  sec  . 

i.e..  when  the  stagnation  zone  is  at  the  target.  At  this  minimum.  \Pt3ieex)  =  0.07  MW , 
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which  corresponds  to  a  propagation  efficiency  of  <  5%.  Increasing  the  slew  rate  further 
brings  the  stagnation  zone  closer  to  the  laser  source,  where  the  intensity  is  lower,  and 
increases  the  effective  wind  velocity  near  the  target.  The  result  is  that  the  power  on  target 
increases  with  slew  rate.  For  a  slew  rate  of  6  =  5  mrad/sec ,  (P^ )  %  0.7  MW .  These 

results  are  consistent  with  the  calculation  in  Sec.  II,  i.e.,  for  a  beam  focused  onto  the 
target  plane,  a  stagnation  zone  close  to  the  target  will  result  in  the  largest  beam  distortion. 

Figure  13  plots  the  laser  fluence  on  target  calculated  over  a  1  sec  dwell  time. 

Panel  labels,  (a),  (b),  (c),  and  (d)  correspond  to  the  different  slew  rates  indicated  in  Fig. 
12.  Panel  (a)  corresponds  to  a  slew  rate  Q  =  —5  mrad/ sec ,  which  gives  the  maximum 
power  on  target  over  the  range  of  slew  rates  considered.  The  beam  is  well-focused  and 
exhibits  practically  no  thermal  blooming  effects.  Panel  (b)  shows  the  beam  profile  when 
the  slew  rate  is  zero,  i.e.,  a  uniform  wind  profile.  In  this  case  the  beam  fluence  profile 
shows  the  characteristic  crescent  shape  associated  with  thermal  blooming  with  the 
centroid  of  the  beam  deflected  opposite  to  the  direction  of  the  wind.  Panel  (c)  shows  the 
beam  profile  on  target  for  the  case  when  the  stagnation  zone  is  at  the  target.  The  beam  is 
severely  defocused  as  a  result  of  thermal  blooming.  The  peak  in  fluence  near  (x  =  0,  y  = 
0)  is  due  to  the  leading  edge  of  the  beam  (t  <  20  msec),  which  is  not  affected  by  the  time 
dependent  thermal  blooming  process.  Panel  (d)  shows  the  beam  profile  when  the 
stagnation  zone  is  1  km  from  the  target.  In  this  case,  the  effective  wind  velocity  near  the 

target  is  in  the  opposite  direction  relative  to  that  of  panel  (b).  The  value  of  (P^ )  is 

similar  to  that  of  panel  (b)  but  with  the  beam  deflected  in  the  opposite  direction. 

Figure  14  shows  the  time  dependence  of  the  beam  intensity  at  the  target.  Figure 
13a  corresponds  to  a  slew  rate  of  6  =  -5  mrad/sec ,  when  the  power  on  target  is 
maximum.  In  this  case  the  beam  is  well-focused  over  the  entire  dwell  time.  The 
transverse  profile  is  mainly  affected  by  turbulence  and  pointing  jitter.  Figure  13b 
corresponds  to  a  slew  rate  of  <9  =  1  mrad/ sec ,  which  results  in  the  minimum  power  on 
target.  For  this  case,  it  is  seen  that  severe  defocusing  occurs  for  times  >  20  msec,  while 
the  first  -10  msec  of  the  beam  is  relatively  unaffected  by  thermal  blooming. 
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Figure  13:  Laser  fluence  profiles  at  the  target  range  for  slew  rates  (a)  6  =  -5  mrad/sec 
(b)  8  =  0 ,  (c)  0  =  1  mrad/sec ,  and  (d)  9  =  5  mrad / sec  .  Parameters  are  the  same  as 
those  of  Fig.  1 1 .  Fluence  is  calculated  over  a  1  sec  dwell  time. 
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Figure  14:  Laser  intensity  in  the  y=0  plane  versus  transverse  coordinate,  x,  and  time.  I,  for 
slew  rates  (a)  6  =  -  5  mrad  /  sec  and  (b)  d  =  1  mrad  /  sec .  Parameters  correspond  to 
those  of  Fig.  1 1. 


VIII.  Conclusions 

In  this  study  we  have  analyzed  the  physical  processes  that  affect  the  propagation 
of  high  energy  laser  (HEL)  beams  and  employed  HELCAP.  a  fully  three-dimensional, 
time-dependent  numerical  simulation  code,  to  determine  the  optimum  laser  wavelength 
and  power  for  HEL  propagation  in  maritime,  desert,  rural,  and  urban  environments.  The 
aerosol  absorption  and  scattering  coefficients  which  characterize  these  environments 
were  generated  using  the  AN  AM  and  MODTRAN  aerosol  models. 

The  theoretical  model  and  numerical  simulations  contain  several  interrelated 
physical  processes  which  affect  HEL  propagation.  These  include:  i)  aerosol  and 
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molecular  scattering  ii)  aerosoi  heating  and  vaporization  iii)  thermal  blooming  due  to 
both  aerosol  and  molecular  absorption  iv)  atmospheric  turbulence  and  v)  laser  beam 
quality.  HELCAP  is  unique  in  that  it  contains  all  of  these  physical  processes  in  a  fully 
time-dependent  and  self-consistent  manner. 

Using  HELCAP.  we  calculated  the  average  power.  ^Pma ) ,  delivered  to  a 

100  cm2  cross  sectional  area  target  at  a  range  of  5  km  over  a  I  second  dwell  time  as  a 
function  of  transmitted  power.  PT .  and  wavelength.  We  considered  three  laser 
wavelengths  corresponding  to  molecular  (water  vapor)  transmission  windows.  1 .045  pm . 

1 .625  pm  and  2.1 41  pm .  and  transmitted  powers  PT  up  to  3  MW.  We  note  that  in 
addition  to  the  power  propagated  to  the  target.  \Pmt, ) ,  the  absorption  efficiency  of  the 

target  should  be  considered  in  evaluating  HEL  lethality.  Target  absorption  efficiency  is 
relatively  insensitive  to  wavelength  in  the  range  considered. 

We  find  that  aerosols  are  of  particular  importance  because  they  result  in  laser 
scattering,  absorption,  and  enhanced  thermal  blooming.  In  the  water  vapor  transmission 
windows,  the  total  absorption  coefficient  driving  thermal  blooming  can  be  due  mainly  to 
aerosols  and  not  to  water  vapor.  In  certain  environments  and  for  sufficiently  high  laser 
power,  scattering  and  absorption  by  aerosols  can  be  reduced  by  vaporization.  Aerosols 
that  consist  of  dust.  soot.  etc.,  cannot  be  vaporized  and  can  significantly  enhance  thermal 
blooming.  We  note  that  moderate  values  of  the  laser  beam  quality  factorM2.  i.e..  values 
less  than  4.  have  little  effect  on  the  propagation  of  HELs  compared  to  other  effects  such 
as.  molecular/aerosol  thermal  blooming  and  turbulence. 

Our  results  show  that  the  average  power  on  target  is  strongly  dependent  on  the 
atmospheric  environment  as  discussed  in  the  Summary.  When  the  propagating  efficiency 
is  not  a  particularly  sensitive  function  of  wavelength,  as  in  the  maritime  and  desert 
environments  for  P1  <  1 .5  MW',  other  issues  such  as  laser  availability  and/or  eye  safe 
wavelengths,  may  become  important  considerations  in  determining  the  optimum 
wavelength  and  powder.  In  this  study  we  have  reported  only  on  wavelengths  in  the  water 
vapor  transmission  windows,  A  =  1 .045 ,  1 .625,  and  2.141  jim .  Other  laser  wavelengths, 
such  as  those  particular  to  solid  state  and  chemical  lasers,  have  also  been  considered.  In 
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general,  it  is  found  that  at  high  laser  powers,  for  which  thermal  blooming  is  a  factor, 
operating  in  the  water  vapor  window  results  in  higher  propagating  efficiencies.  We  find 
that  propagation  efficiencies  approaching  60%  can  be  achieved  in  the  absence  of 
stagnation  zones. 

Experiments  are  presently  underway  at  the  Naval  Research  Laboratory 
to  study  and  characterize  aerosol-induced  thermal  blooming  and  scattering  of  HEL 
beams.  In  these  experiments  a  -  1  kW,  CW.  -  1pm  fiber  laser  will  interact  with  both 
water  based  and  dust  like  aerosols.  Earlier  experiments  performed  at  United  Aircraft 
Research  Laboratories  using  and  a  low  power  15W.  10.6pm  laser  beam  demonstrated 
enhanced  thermal  blooming  in  the  presence  of  carbon  aerosols  [7]. 

Stagnation  zones  can  significantly  degrade  the  propagation  efficiency  of  a  high 
energy  laser  (HEL)  through  the  atmosphere.  We  have  analyzed  the  propagation  of  a 
focused  HEL  beam  through  a  model  atmosphere  containing  a  stagnation  zone.  Our 
analysis  shows  that  the  amount  of  laser  degradation  caused  by  a  stagnation  zone  is  highly 
dependent  on  the  focusing  geometry  of  the  laser  beam  and  the  location  of  the  stagnation 
zone  along  the  propagation  path.  For  a  collimated  laser  beam,  a  stagnation  zone  near  the 
transmitter  causes  the  largest  distortion  of  the  laser  beam  on  the  target.  For  a  laser  that  is 
focused  onto  the  target,  a  stagnation  zone  near  the  target  causes  the  largest  laser 
distortion  and  the  fastest  decrease  in  laser  intensity.  This  finding  is  in  agreement  with  the 
experimental  results  of  Ref.  [16]. 

We  used  the  HELCAP  to  model  the  propagation  of  a  megawatt-class  HF.L  beam 
through  a  realistic  maritime  atmosphere  containing  a  stationary  stagnation  zone.  The 
HEL  beam  is  focused  onto  a  remote  target.  The  stagnation  zone  is  created  by  slewing  the 
HEL  beam  in  the  direction  of  the  wind  and  location  of  the  stagnation  zone  is  varied  by 
changing  the  slew  rate.  Consistent  with  our  analysis,  the  simulations  show  that  the  power 
on  target  is  minimized  when  the  slew  rate  is  such  that  the  stagnation  zone  is  located  near 
the  target  plane.  In  this  case,  propagation  efficiency  can  be  reduced  by  an  order  of 
magnitude  relative  to  an  unslewed  beam. 
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